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Abstract. The Orion Nebula is one of the most frequently observed nearby (< 1 kilo- 
parsec) star forming regions and, consequently, the subject of a large bibUography 
of observations and interpretation. The summary in this chapter is bounded spatially 
by the blister Hll region, with sources beyond the central nebula that are part of the 
same dynamical clustering covered in other chapters in this book. Herein are discussed 
panchromatic observations of the massive OB stars, the general T Tauri population, the 
sub- stellar sources and variable stars within the Orion Nebula. First, a brief history of 
400 years of observation of the Nebula is presented. As this history is marked clearly 
by revelations provided in each age of new technology, recent ultra-deep X-ray surveys 
and high resolution multi-epoch monitoring of massive binary systems and radio stars 
receive special attention in this review. Topics discussed include the kinematics, mul- 
tiplicity, mass distribution, rotation, and circumstellar characteristics of the pre-main 
sequence population. Also treated in depth are historical and current constraints on the 
distance to the Orion Nebula Cluster; a long standing 10-20% uncertainty has only re- 
cently begun to converge on a value near ^ 400 parsecs. Complementing the current 
review of the stellar population is a companion chapter reviewing the molecular cloud, 
ionized Hll region and the youngest protostellar sources. 



Introduction 

An interesting hypothesis drawn from our knowledge about the Orion Nebula is that 

50, 000 years ago it was invisible to the naked eye. The ionizing photons of the massive 
O and B type stars, whose projected arrangement yield the namesake Trapezium, had 
not yet burned away the layers of natal molecular gas out of which they had formed. 
While bright blue stars were visible along the Sword of Orion, having formed continu- 
ally over the previous few milhon years, there were, on a scale perhaps much grander 
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Figure 1 . Color composite image of the northern Orion Molecular Cloud. This 
ground based image is oriented with equatorial North up and East to the left and 
encompasses the central area of the young star distribution shown in Figure|2] Image 
courtesy Robert Gendler. 
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than the present, thousands of smaller stars hidden from view. In the relatively short in- 
terim a blister HII region created by a newborn 40, OOOK O type star expanded into the 
molecular cloud, uncovering a large portion of this embedded star clustering. Nonethe- 
less, star formation continues vigorously in the remaining molecular cloud today. 

Because of the richness of this star clustering (A^^ > 2000) and its relative prox- 
imity (~ 400pc), the Orion Nebula is easily the most frequently observed nearby 
(< 1 kpc) star forming region, providing an large bibliography of observations and 
interpretation^ Moreover, the properties of the Orion Nebula stars, e.g., their masses, 
evolutionary status, spatial and velocity distributions, outflows and circumstellar disk 
properties, all provide critical tests for theories of molecular cloud evolution and star 
formation. 

Region Overview 

Inspecting a magnificent modern large scale visual image of the Orion Nebula (Fig- 
ure [T]) reveals the major physical features of this region. From North to South there are 
a series of bright emission nebulae, interspersed with dark bands and small clusterings 
of bright stars. The Orion Nebula HII regiorj^is central to this image and appears to 
expand to the southwest from an apex at the location of the O and B stars. These appar- 
ent alternating nebulae and clusterings have led to a system of names or designations 
with boundaries that deserve some explanation. 



Table 1 . Subregion names in the northern Orion A molecular cloud. Papers in- 
cluded : [PEenago|ll954a| ) ;|Blaauw| ( fT964ll; [Morgan & Loden|(|T966| M&L); |Waiker 



( fT969] l; |WarTCn & Hesser| ( fr978| W&H)I lGomez & Lada| ( |1998| G&L). Note that the 
G&L98 subclusters F and G, which are outside the immediate field of study, appear 
to correspond to the LDN 1641N and NGC 1999 star forming regions, respectively, 
although the G&L98 Ha subclusters are shifted ^ 10' to the West. These regions 
are discussed further the chapter by Allen & Davis. 



Adopted Name 


Parenago 


Blaauw 


M&L 


Walker 


W&H 


G&L 




1954 


1964 


1966 


1969 


1978 


1998 


Upper Sword 


I 


Orion Ic 


Upper Sword 


Group 1 


CI 


A? 


NGC 1977 


II 


Orion Ic 




Group 2 


C2 


A 


OMC 2/3 


III 


Orion Ic 




Group 3 


C3 


B 


ONC 


IV 


Orion Id 




Group 4 


D/Dl 


CD 


iOri 


V 


Orion Ic 


tOri 


Group 5 


C4 


E 



As listed in Table [T] the stars along the Sword of Orion have traditionally been 
segregated into 4 or 5 regions; as shown in the tabulation the "names" for these re- 
gions have changed over time although no new divisions have been made since the 
1950s. Subsequently, the numeral ordering (I,II,...) from Parenago (1954) is adopted 



The introductory sketch of the recent history for the Orion Nebula is but one hypothesis taken from the 
breadth of observational and theoretical studies of the Nebula. A more complete summary of such models 
for this history are presented in a companion chapter (O'Dell et al). 

^Additional common catalog entries for this region include Messier 42, NGC 1976. 
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Figure 2. Young star distribution along the Orion Molecular Cloud. The variable 
infrared stars identified by Carpenter et al.| ( |2001 )(yellow circles) are compared to 
the OB members of the Orion Ic (blue diamonds) and Id (red squares) ( |Brown|1996|l. 
The reverse grayscale image is the MSX band A (8 micron) image from 'Kraemerj 
et al. I ( [2(303 1. Labels (I. ,11.,...) mark the five subregions of Parenago (Table [T]i and 



their more descriptive names. The equatorial coordinate axis (decimal degrees) are 
in equinox J2000. This chapter focuses on the stars forming in region IV, the Orion 
Nebula. 



and paired with more descriptive names, e.g.. Region II corresponds to the NGC 1977 
HII region. 

While at optical wavelengths these features seem apparent (and perhaps some such 
as NGC 1977 are significant) when one views the Sword of Orion at near-IR wave- 
lengths (2 fim), which are much less sensitive to variations in extinction, no boundaries 
are apparent between these regions. Figure [2] presents a map similar in extent to the 
previous optical image but where the Nebula is traced by mid-IR (8 //m) emission as 
observed by the Midcourse Space Experiment (MSX) ( Kraemer et al.|2003 1. Infrared 
variable stars ([Car penter et al.|200T I extend across all these groups without clearly de- 
marcating any of them except for the strong concentration of sources at the apex of the 
Orion Nebula. Another set of frequently used designations concern the OB stars of the 
Sword; [Blaauw [ ([1964|) segregated the Sword OB stars in the Ic and Id associations, 
ordered in part by apparent youth, while Warren & Hesser ( 1978 ) expanded the Ic asso- 
ciation on the sky but segregated the Sword OB stars into subregions (CI, C2, C3, C4), 
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which are referred to in this review by the collective designation Ic*. While the Orion 
Id OB "members" coincide with the stellar density maximum in Figure |2j the spatial 
division of Ic* and Id members appears rather arbitrary. 

Another naming system deserves clarification. The terms "Trapezium" cluster 
([Trumpler'1931), which refers to stars immediately surrounding the asterism that is the 
arrangement of 4 bright OB stars in the center of the Nebula, and "Orion Nebula Cluster 



(ONC)," which dates to Haro (19531, suggest perhaps that these are separate entities. 
This survey of the literature does not reveal any physical reason to suppose that the 
Trapezium stars represent anything distinct about star formation in the Nebula beyond 



their mass. While the study by Hillenbrand ( 1997 1 arbitrarily divided the region into 



three "radial zones:" Trapezium Cluster (r< 0.3 parsec); ONC (r< 3 parsec); Orion Ic 



association, only slight age gradients between them were found (see also|Ramfrez et al. 



2004 1. It is therefore a secure inference that the entire region is a single contiguous star 



forming event that requires complete description. 
Topical Scope of the Review 

This chapter summarizes current knowledge regarding those stars which sit within and 
surround the Orion Nebula HII region. A review of the cold molecular cloud, the hot 
HII region and the characteristics of the very youngest stars is presented in a subsequent 
chapter (O'Dell et al.). Further outlying regions in Figure [2] e.g. to the North (OMC 
2-3; NGC 1977, Upper Sword) and South (LDN 1641) of the Nebula are discussed in 
other chapters in this volume (i.e., Peterson & Megeath and Allen & Davis, respec- 
tively). The focus here is on those surveys that provide consti^aints on the physical 
properties of the revealed high mass, T Tauri and substellar objects. 

The review is organized as follows. First, a brief history of important or previ- 
ously broad reviews of research on the Orion Nebula is presented (Sect. [T]]l, followed 
by an overview of distance determinations to the Nebula in the past 70 years (Sect.[2!]). 
Individual sections are reserved for reviewing the T Tauri members with a focus on 
summarizing the many broadband CCD and spectroscopic studies that have occurred 
during the past two decades (Sect.|3]]l, the O and B type members of the Nebula with de- 
tailed reviews of each of the Trapezium stars (Sect.|4]]l, and the variable stars (Sect. [5^. 



1. History of Study of the Orion Nebula 



During the fifty years after the development of the telescope in 1608 the nebular nature 
and stellar content of Orion Nebula were independently discovered by a handful of 
observers. The observing logs of Nicolas-Claude Fabri de Peiresc (1610) and of Johann 



Baptist Cysat (at the latest 1618) ( |Wolf||1854] |Holden|[T882] ) represent the earhest 
written records that the Sword of Orion contained a "fog" or "milky nebulosity," which 
here borrows the words used by W. Herschel (1802) to describe the region. The first 
hand drawn charts of the region include those of Galileo (1617), who did not distinguish 
the Nebula, of Giovanni Battista Hodierna (drawn sometime before 1654), who did, and 
the more famous 1656 drawing by Huyghens, which is the most widely known. His 
accurate rendition of the central nebula surrounding the Trapezium provides an origin 
for the term "Huyghenian region" (Figure |3]l. 

A nearly complete journal of 273 years of telescope aided visual observation of 
the Nebula is provided in Edward S. Holden's The Monograph of the Central Parts of 
the Nebula of Orion (1882). Holden's monograph includes the reproduction of dozens 



Figure 3. The 1656 sketch of the inner Orion Nebula by Huyghens. The figure 
has been re-oriented from its publication form such that North is up and East is to 
the left. Image is reproduced from |Holden| ( | 1 882) l . 



of hand drawn sketches as well as observing logs each in their original language. Vari- 
ation in the reproduction of the Nebula is remarkable. The first photographic plates of 
Orion made by Henry Draper between 1880 and 1882 were included as an addendum 
to Holden's work as well as a discussion of the processes of obtaining these images. 
Figure [4]is a reproduction of that image and is captioned with Holden's description. 

The next 100 years of photographic observation of the Nebula included quanti- 
tative studies of its variable stars, the discovery of a cluster of faint stars in the Nebula's 
core, very broad censuses and an expansion in the role of the Nebula as a testbed for 
new observing techniques. Numerous variable star studies were performed using the 
Harvard Plates by H. Leavitt (published by Pickering) or confirmed by M. Applegate 
(published by Shapley). In the 1930s, deep red photographic plates revealed that the 
Nebula contained a substantial cluster of fainter stars in addition to the brightest mem- 
bers JTmmpeFTWT;'Baad£^&^^^wsE][^^ The broad surveys of Brun (1935l 
and [Parenago, (1954a) provided excellent photographic updates to the ,Bond| (1867 1 
visual census of stars in the Nebula. 

Parenago, a Russian astronomer, published a major, lengthy analysis of the region 
in 1954. A translation of sections of this publication was undertaken for this review and 
it indicates that he relied on prior and concurrent work of female Russian astronomers 
(e.g., Barkahatova, Uranova, Kirillova) some of whose names do not reappear in the 
literature outside of his book. His analysis extended to the topics of parallax, proper 
motions, astrometry and photometry, covering all of the Sword of Orion. Most impor- 
tant he found a clear evidence for a "cloud" of members lying above the main sequence. 
Perhaps because of the Cold War and the lack of translations from his work from Rus- 
sian to other languages, his work is extremely poorly cited in the literature. The dis- 
paraging of his work by Walker (1956 1 was not unnoticed by the author as revealed 
in the posthumous publication, Parenago & Sharov ( 1961| ). Nevertheless, this work 
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Figure 4. Draper's 1882 photographic image of the Orion Nebula. North is up and 
East is to the left. An approximately 2 hour visible photograph obtained by H. Draper 
with the 11 inch Clark telescope. An excerpt from Edward S. Holden's "Monograph 
of the Central Parts of the Nebula of Orion" (Washington 1882) p. 226: "The first 
photograph of the nebula of Orion was made by Dr Henry Draper in September, 
1880, and the unavoidable delay which has occurred in printing the present memoir 
enables me to include an account of the astonishing results which he has attained. A 
wood-cut which I had prepared from his first photograph was found to be so unsat- 
isfactory that Dr. Draper most generously offered to supply the necessary photolito- 
graphic reproductions of his last negative (taken March 14, 1882) to accompany the 
brief account I had prepared. The full page photolitograph is here given..." 



was not considered seriously in most subsequent studies, not appearing, for example, 
in the otherwise meticulous study by Goudis (see below). Even today this significant 
work has garnered a mere 59 citation ^ H is valuable data tables were converted to 
machine readable formats by Malkov ( 1992 ) and ingested into an electronic format in 
1997 ( |Parenago|1997l ). 

In 1982, approximately one hundred years after the Holden monograph and the 
first photographic images by Draper, two useful summarizing publications appeared. 
C. Goudis's The Orion Complex: A Case Study of Interstellar Matter, focused on the 
structure and nature of the Nebula itself, details previous approaches to studying the 
Nebula and its content, including infrared, radio and spectroscopy, and provides use- 
ful tables of past observations and results. Second, a conference was held on the Orion 
Nebula and honoring Henry Draper ( Glassgold et al. 1 1 982 1 ; the conference proceedings 



^Derived via the NASA Astrophysical Data Service circa May 2008 by merging the results of citations to 
Parenago 1 1954a^ and Parenago 1 1954b i. 
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include 32 articles on all current aspects of study of the Nebula and records of the par- 
ticipants discussion about each contribution. In addition there are a number of articles 
that provide a history, more detailed than that of Holden, about Draper's photographic 
work obtaining these images of the Nebula as well as his subsequent scientific legacy 
after his death in late 1882. 

Reviews of the literature over the subsequent ~ 20 years of CCD observations of 
the Nebula include [GeiTzel & Stutzki] ( fT989] l and |O'Dell| ( [200T] ). The Nebula was also 



included in a review that encompassed all of the Orion star forming region by Brand & 



Wouterloot ( 1991 1. Its recent study has also been the focus of a book, ,0' Dell, ( ,2003 1. 



2. Distance 



Minkowski ( 19461 opens his re-analysis of the Trumpler] (|1931') derivation of the dis 



tance to Orion with a sentence still applicable today, "All published values of the dis- 
tance of the Orion Nebula are open to some criticism." At the time, the range in quoted 
distances was a factor of 10. While the spread in acceptable values has decreased over 
time, uncertainty in the distance to the stars in or near the Orion Nebula at the 10-15% 
level remains today. 

The difficulty in estimating distances is due in part to the complex geometry and 
kinematics of the region as a whole, and in part to characteristics of the youthful mem- 
ber stars themselves. Several stellar subgroups were identified by Blaauw|(|1964|), each 



covering several degrees on the sky. These groups, which appear to have different ages, 
overlap along the line of sight with a total depth of more than 100 pc. This renders 
membership boundaries and hence distances to the individual subgroups (as well as 
ages) difficult to distinguish. The molecular cloud containing the Orion Nebula Cluster 
(ONC) is behind most of the optically visible early type stars in the larger association. 
Depending on the sub-group within the association, the O and some of the B stars are 
slightly evolved from the ZAMS while the A and/or later stars may be still in the pre- 
ZAMS phase. Reddening is spatially variable and significant, especially towards the 
Nebula. Further, the vast majority of nebula stars are photometrically variable and have 
other signatures of circumstellar activity in their photometry and spectra. 

Table |2] contains relevant distance determinations to the Orion Nebula Cluster and 
records, where possible, the method used to derive the distance and if error estimates 
were documented. Authors were found to have often included some mix of Id, Ic*, and 
Ic stars in their distance estimates; occasionally authors provided estimates derived for 
the outer parts of the Id association but applied to the inner region. Therefore Table [2] 
includes all Ic and Id distance estimates; because it is well established that Ic members 
experience less line of sight extinction than their Id counterparts it is secure to infer that 
Ic distances provide lower limits to the distance to the cluster. 

Four basic methods provide most of the distance estimates to the Orion region: 
zero-age main sequence fitting or similar stellar evolutionary status methods that pro- 
vide distance moduli, kinematic methods that assume a specific dynamical model for 
the cluster, parallax estimates, and reddening analyses. Most authors have used B stars 
as distance probes since they are bright (enabling good data to be obtained) and close 
to the ZAMS (enabling distance modulus determination). Below, each distance estima- 
tion method is discussed in turn, followed by a synopsis of the best current constraints 
on the ONC distance, and finally a survey of commonly cited references for the ONC 
distance that, in fact, do not contain actual analyses of the ONC distance. 
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Table 2. Summary of published distances to the Orion Id association. Distances 
are segregated as corresponding to the Ic or Id associations with the notation Ic* 
used to indicate that part of Ic right around Id and to distinguish it from the much 
larger Ic as defined by WH78. Unfortunately the notation Ic* was also used by many 
authors before the Ia,b,c,d subgroups were defined to indicate the entire "Sword" 
region; thus some ambiguities may remain. 



Yeai- of 
Pub. 



Region Distance Error 
Desig. Modulus 



Stellar 
Types 



Nuinber 
of Stars 



|Pickerin^ 

Kaptevn 
nciierinu 
^'annei^oel^ 
rrumpier 
iviinicowsici 
^iviari^owitZj 
janarpressj 



1917 
191" 
191 

1929 
1931 
1946 
1949 
1952,4 



/ 



Id+Ic* 

Ic 
Id+Ic* 

all 

Id 

Id 

Id 
Id+Ic 



11.5 
6.34 
8.5 
7.6 
8.5 
7.38 
8.58 
8.5 



0.35 
0.30 



B3 

B 
B3 

B 

09-A2 
07-Bl 
B1-B3 
B Stars 



ZAMS 

PM 
ZAMS 
ZAMS 
ZAMS 
ZAMS 
ZAMS 
ZAMS 



Pg 
Pg 
pg 
pg 
pg 
pg 
pg 
pc 



17 

3 
17 
190 



I Parenagol 



UmillSUIl &. Hiltnerl 



itrantl 
Stiarpiess 

Boriiman & Blaauw 
Jonuson 

iviorsian & Loden 

Lesh 
I Walker! 
[Peristbri 
[P^ristorl'et al.| 



1952 


Id 


8.6 


0.30 


B Stars 


ZAMS 




1 


1954 


Id 


8.0 






ZAMS 


pg 


7 


1956 


Id 


8.0 




B8-A0 


ZAMS 


7 


1958 


Id 


8.6 




06-K2? 


PM/RV 


plate 


20 


1962 


Id 


8.2 




B Stars 




UBV, H-, 


180 


1964 


Id 


8.33 


0.11 




ZAMS 


7-filter 


5 


1965 


Id 


7.9 






PM/RV 




21 


1966 


Ic* 


8.1 






ZAMS 




7 


1968 


Ic* 


8.5 


0.1 




ZAMS 




14 


1969 


Id+Ic* 


8.37 


0.05 


B2-B9 


ZAMS 


UBV 


51 


1973 


Id 


7.8 


0.15 


B stars 


ZAMS 


BV 


15 


1975 


Id 


8.1 


0.13 


B stars 


ZAMS 


BV 


7 


1975 


Id 


7.71 


0.21 


B stars 


ZAMS 


VI 


7 



. 1975 

IWarren & Hesserl 1977 



7.98 
8.42 



0.12 
0.53 



B stars 
B stars 



V;SpT ? 
UBV; ubvy; Hfj 6 




1977 
1981 
1981 
1981 
1982 



Ic* 
Id+Ic* 
K-L region 
Id+Ic 
Ic 



8.16 
8.20 
8.41 
8.0 
7.87 



0.49 
0.15 



ZAMS? 



0.40 H20masers stat. parallax 



0.5 
0.09 



ZAMS 
ZAMS 



UBV; ubvy; Hf^ 
UBV? 
VLBI pm + rv 
BV 
Mv;H„ 



1982 
1990 



Id+Ic* 
Ic 



8.19 
7.7 



0.10 
0.50 



i stars 
i Stars 



ZAMS 
ZAMS 



Mv;Hp 

CO-H-y 



[Brown et al.| 



1990 
1994 



8.2 
8.0 



0.03 
0.49 



i Stars ZAMS 
i Stars ZAMS 



VBLUW 



2 
34 



Brown et al 
de Zeeuw ei 
.Srassun e[ 



ei afl 
alf*. 
H6i^dridgz^al.| 



1- 

[Kraus et al.| 



[Hirota et aTl 
"SandsironT^ al.l 
ivienien ei ai. ( 
iviayne 6l iNaylor 



1994 


cloud 


8.1 


0.48 


B Stars 


Red. 


VBLUW 




1998 


Ic 


8.32 


0.17 


B Stars 


trig, parallax 


Hipparcos 


34 


1999 


Ic 


8.52 


0.25 


B Stars 


trig, parallax 


Hipparcos 


34 


2004 


Ic* 


7.96 


0.10 


ec. binary 


Radius 




1 


2005 


Ib+Ic 


8.23 


0.08 


B5-F Stars 


trig, parallax 


Hipparcos 


121 


2005 


Ib+Ic 


7.97 


0.10 


B5-F Stars 


ZAMS 


BV/Hipparcos 


111 


2005 


Id+Ic 


8.34 


0.32 


Stars 


Red. 


CO/Hipparcos 


7 


2006a 


Ic* 


8.19 


0.30 


ec. binary 


Radius 




1 


2007a 


Id 


8.22 


0.16 


G6-M2 


R sin{i) 


Various 


74 


2007a 


Id 


7.97 


0.17 


G6-M2 


R sin{i) 


Various 


34 


2007 


Id 


8.19 


0.06 





dyn. parallax 


binary orbit 


1 


2007 


Id 


7.94 


0.06 





dyn. parallax 


binary orbit 


1 


2007 


K-L region 


8.20 


0.09 


H2O masers 


trig, parallax 


VERA 


1 


2007 


Id 


7.95 


0.13 


radio star 


trig, parallax 


VLBA 


1 


2007 


Id 


8.08 


0.03 


radio stars 


trig, parallax 


VLBA 


4 


2008 


Id 


7.96 


0.06 


Bl-AO 


MS Fitting 


VI 


~20 
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2.1. Distances Based on Stellar Evolution 



Assumption of Zero-Age Main Sequence Stars Spectroscopic parallax distance esti- 
mates involve comparison of (de-reddened) apparent magnitudes with absolute mag- 
nitudes, which are assumed based on stellar evolutionary status, to derive a distance 
modulus. Several early works (see Table [2]l reported distances ranging from 185 to 
2000 pc (falling off the range of distance versus publication year shown in Figure |6]). In 
the first modern set of distance estimates to the Orion Nebula region, Trumpler] derived 
a distance modulus (DM) of 8.5 magnitudes by comparing absolute magnitudes (from 
|Trumpler| |1930) along the main sequence with the de-reddened apparent magnitudes 
of B stars in the cluster]^ His method would be repeated by numerous authors in later 
studies, |Minkowski re-calculated the total absorption towards the three brightest stars 
in the Trapezium, used a different absolute magnitude relation (that of Wilson||1940 1, 
and derived a much smaller distance modulus, 7.38 magnitudes. These first two modem 
distance estimates are the extremes in distance found in Table |2] 

In later work, Markowitz ( 1949 1 similarly studied early B stars, used yet another 
absolute magnitude relation (that of Blaauw 1946), and found a distance modulus of 
8.57. Sharpless (1952i performed the most extensive survey yet of early type stars 
throughout the entire Orion region with special focus on stars within a few degrees of 
the Trapezium. He found a distance modulus of 8.5 for the ensemble and a slightly fur- 
ther 8.6 for stars near the Nebula. Using the "Q-method" of photometric dereddening. 



rather than spectroscopy, and the Johnson & Morgan ( 1953 1 main sequence, Sharpless 



( 1954 1 reported the same distance modulus of 8.5. Parenago ( 1954a I derived a distance 
modulus of 8.0 from his vast photographic catalog; he used A stars rather than the B 
stars typical of other authors. 

Johnson & Hiltner'(|1956|l recognized that some luminosity evolution away from a 



zero-age main sequence will occur between clusters of different ages and that this may 



be the case for the B stars in the ONC. Using a re-calibration of the Johnson & Morgan 



( 1953 1 main sequence, they calculated a distance modulus of 8.0 by de -reddening the 



data of Sharpless ( 1954 1 for stars below the assumed upper MS turnoff. However, their 
assertion that A stars are on the main sequence rather than in the pre-main sequence 
phase of stellar evolution is likely not correct for the young Orion Nebula Cluster, lead- 
ing to an underestimate of the distan ce. |Sharpless, ( ,1962) used photoelectric data to 
revisit and revise the Sharpless (1952 1954 1 distance downward to 8.2 magnitudes by 
considering the evidence for stellar evolution. Additional applications of traditional 
BV or U BV zero-age main sequence fitting include IBorgman & Blaauw|(|1964[),|Mor- 



gan & Loden| (1966), Lesh (1968), Walker (1969) and Pension (1973) - revised by 
Pension et al. |1975 ). Herna ndez et al.| ( |2005) ) used Hipparcos photometry and the main 



sequence of [Cox, ((2000). These authors contended with reddening in different ways and 
selected to varying extents samples free of binaries or variable stars. 

Warren & Hesser| ( |1978 1 first applied narrow band Stromgren photometry of B 



stars to the Orion distance problem, deriving a frequently cited distance to the Ori 
Id region of 435 pc, as well as distances to the other subgroups and newly defined 



sub-divisions of the subgroups. Anthony-Twarog ( 1982 1 revised the Warren & Hesser 
distance using the same data but a different H/? calibration and different combinations 



*The result of a second method based on the angular diameter of the cluster was also reported. He assumed 
the Orion cluster had a size typical of other open clusters. 
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of subgroups. In general, th e |Anthony-Twarog| distance estimates are 40-80 pc closer 
those of Warren & Hesser [Wolff ( 19901) used H7 and the Balmer discontinuity to 



determine Tg//, surface gravity and the absolute bolometric magnitude of B stars in 
Orion, also deriving distance estimates to all four subgroups in Orion. 

Brovyn et al.| ([1994') used VLU BW photometry and interpolated grids of Ku- 



rucz models to derive stellar parameters from which ZAMS fitting techniques could 
be applied. Distances to each OB association sub-group were derived and these re- 



sults display a systematic 0.3 mag shift (closer) than those derived by Warren & Hesser 



( 1978J). Mo re recently, |Mayne & Nay lor ( 2008,) use d the photometry and effective tem- 
peratures of Hillenbrand \ 1997 ), the Mathis ( 1990 1 extinction law, and Geneva-Bessell 
isochrones to derive 391 pc; this is closer than previous estimates using main sequence 
fitting techniques but consistent with contemporaneous distance estimates using other 
techniques. 

In summary, distance estimates derived using assumed constraints on the evolu- 
tionary status of the Orion OB stars are widely scattered. The primary uncertainties 
arise from sample selection, reddening corrections, assumptions about the evolutionary 
state of the early type stars, and the adopted main sequence which can vary by several 
tenths of a magnitude between authors. 



Stellar Rotation A characteristic of young stars is their relatively rapid rotation, which 
is measured using time-series observations that track the periodicity of cool spots on a 
star's surface or through spectroscopic measurement of the velocity broadening of ab- 
sorption line features. Coupling these two observations yields a distance independent 
measure of a star's radius convolved with the inclination of the star's rotational axis on 
the sky. Comparing the radii derived from the cluster stars' rotational properties to that 
derived from placing the stars on the HR diagram, one can derive the distance to the 
cluster if one assumes that as an ensemble the stars in a cluster have randomly oriented 
rotation axes the sky. This technique was first developed by Hendry et al.l (|1993 1 and 
has been applied to the Pleiades (O'Dell et al. |1994[ ), and Taurus (Preibisch & Smith 
|I997) ). |Jeffries| ( |2007aj ) applied this method to the ONC, using the large database of 
periodic stars with measured rotational properties. The fact that the canonical distance 
(480 pc) to the ONC did not yield a randomized sin{i) distribution was first shown by 
Rhode et^L] {2001 ), who did not, however, estimate the amount that the cluster's dis 



tance would have to decrease. Jeffries cataloged 74 young stars in the Nebula having all 
of the requisite observations and used them to derive a distance of 440 ib 32 pc. After 
showing that the accreting stars appear to have systematically underestimated luminosi- 
ties (and thus a biased sin{i) distribution), Jeffries derived a distance of 392 it 32 pc 
using a subset of 32 non-accreting young stars in the Nebula. 



2.2. Distances Based on Kinematics 

Proper Motions Kinematic distance estimates involve radial velocity and proper mo- 
tion data combined with kinematic assumptions. In the simplest model of random mo- 
tions, the distance is directly proportional to the ratio of the radial velocity and proper 
motion dispersions. More complex models such as expansion, contraction, or rotation 
can be employed as well. [Strand] ( 1958 1 derived the first distance to the ONC using 
this method. He combined proper motion data with only a few radial velocities and 



used an expansion model to estimate a distance of 520 pc to the ONC region. Johnson 



( 1965 ) presented new radial velocities and used the Strand ( 1958| proper motions of the 
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same stars to derive a smaller distance of 380 pc, assuming random motions. Finally, 
Walker ( 1983[ ) presented new radial velocity data and computed both a radial velocity 



dispersion and a proper motion dispersion for the same stars from the data of Parenago 
([1954a I but did not carry the analysis through to a distance estimate. 

The most frequently cited distance to the Orion Nebula Cluster comes not from 
study of stellar motions, but from proper motions and radial velocities of H2O masers in 
the Kleinmann-Low (K-L) nebula by |Genzel et al. ( 1981 1. The K-L nebula is embedded 
in the Orion Molecular Cloud behind the Orion Nebula, though it is thought to be within 
1 pc of the front edge of the cloud (e.g. Zuckerman 1973 1. Maser velocities were 
compared with a kinematic expansion model for the outflow to derive a distance of 
480 ± 80 pc. That expansion model has undergone subsequent changes in its inferred 
orientation on the sky ( [Greenhill et al. 1998, |2004b^ but the impact of these model 
changes on that distance estimate have not been quantified. 

The primary uncertainties in these kinematic methods lie with sample selection, 
with assumptions of the kinematic models, e.g. random motion versus expansion, con- 
traction, or rotation, and with their use over stellar groupings large enough to be consid- 
ered unbound associations rather than bound clusters. The scatter in distances derived 
from kinematic methods is comparable to that in distances derived from zero-age main 
sequence fitting. 



Double-line Eclipsing Binary Systems Further distance estimates can be derived from 
kinematic analysis of double-line eclipsing binary systems, which provide empirically 
constrained values of the stellar radii. In the Orion Nebula plus Ic* association there are 
currently two such systems that have refined results. Stassun et al. ( 2004b I derived a 
distance of 419 ± 21 pc (or 390 pc adopting a more conventional value for the bolomet- 
ric magnitude of the Sun) for VI 174 Ori, an M-type pre-main sequence solar analog 
system. Partly based on this distance and partly based on age arguments, these authors 
consider this star a member of Orion Ic*. A similar, remarkable analysis of the brown 
dwarf - brown dwarf eclipsing binary 2MASS J05352 184-0546085 yielded a distance 
of 435 ± 55 parsecs ( [Stassun et al. |2006a ); both stai^s are projected against the south- 
ern reaches of the Nebula. One systematic that is not constrained by the dynamics of 
these systems is the line of sight extinction; in these cases increasing the inferred line of 
sight extinction acts to move the star's inferred distances to smaller values. Additional 
eclipsing binaries ( Irwin et al.|2007 Cargile et al.|2008 1 will, eventually, lead to further 
fundamental distance constraints. 



2.3. Direct Parallax Determinations 

Hipparcos Hipparcos trigonometric parallax distances to the individual Orion sub- 
groups were first provided by Brown et al. (1998; unpublished preprint; see also dis- 
cussion in Brown et al.|1999 ) who found that Orion la is 50-100 pc in front of the other 
associations (Orion lb, Ic), consistent with the spectroscopic parallax analyses of Brown] 
etaL] ( [T994| ) and |Warren & Hesser| ( |T978 ]). A distance of 462 ± 36 pc (DM = 8.32) was 



quoted for the Ic group, and cited as preliminary. The difficulties in interpreting Hip- 
parcos parallax data in the Orion region of the sky include 1) mostly radial motion of 
both members and field stars due to location towards the solar antapex, which causes 2) 
significant membership biases, while 3) Orion is located close to the upper limit (500 
pc) of Hipparcos sensitivities. 
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The Brown et al. (1998) results were revised by de Zeeuw et al. ( 19991 who also 
acknowledged the astrophysical difficulties of Orion and found a distance of 506ib37 pc 
(DM = 8.52) to the Ic group, using the same stars as Brown et al. (1998). Reasons for 
the different results for the Ia,b,c groups are not explained. Brown et al. ( 1999 1 restated 
the|de Zeeuw et al. results but also noted them as preliminary. 

Finally, Hernandez et al. ( 2005] ) recalculated parallax distances to Orion subgroups 
using a revised B star membership list selected according to kinematic and color cri- 
teria and partially revised sub-group designations. These authors find a distance of 
443 lb 16 pc to the combined Orion lb and Ic regions, in agreement with Brown et al. 
and de Zeeuw et al. within errors. However, the distances to these two subgroups are 
derived together rather than independently and so the implications for the distance to 
Orion Id is unclear. 

The primary uncertainty with existing Hipparcos parallax estimates of the Orion 
distance is its limited precision at such a large distance combined with the astrophysical 
circumstances regarding Orion kinematics. Future missions should redress the first of 
these issues and possibly overcome the second. 



Interferometric Observations of Radio Sources Very long baseline radio interferome- 
try provides the astrometric precision necessary to measure and separate the combined 
proper motion and parallax reflex motion for compact objects at a distance of 500 pc. 
Recent results have utilized the fact that the strong magnetic fields of pre-main sequence 



objects cause them to be excellent radio targets for parallax determinations ( jSandstrom 



et al.|2007 Menten et al.|2007 1, while another group has derived the annual parallax of 



water masers in the Kleinmann-Low nebula ( [Hirota et al.|2007 



Four radio stars in the Orion Nebula have been used to derive distance estimates: 
GMR A, F, G and 12 (GMR: |Garay et al.|[T987l l. [Sandstrom eTal] ( [2007j ) used the 
stai" GMR-A, observed it with the Very Long Baseline Array (VLBA) during the 2003- 
2004 epoch and found a parallax ofvr = 2.57 it 0.15 mas or a distance of 389^21 P'^ 
(Figure [5]). GMR-A is optically obscured by the molecular cloud associated with the 
Nebula, and should, thus, provide an upper limit on the distance to the Nebula. During 
the 2006-2007 epoch 'Menten et al.'(2007) used the VLBA to measure the trigonomet- 
ric parallax of GMR-A and three additionial variable non-thermal radio sources. They 
found a parallax of vr = 2.390 ± 0.104 (418.4 ± 18.2 pc) for GMR-A, and a joint 
solution for all four sources of vr = 2.415 ± 0.040 (414.0 ib 6.8pc). The precision 
of the source positions for such observations are affected by the time variability of the 
targets' flux and the fact that these sources are sometimes resolved (i.e., not concur- 
rently point sources, e.g., GMR-A). In addition these two VLBA studies used different 
calibration techniques. 

Hirota et al. ( |2007 1 observed water masers in the Orion K-L region during the 



2004-2006 period using the VLBI Exploration of Radio Astrometry (VERA) system 
in Japan. Filtering the observed set of water masers based on signal-to-noise and on 
vlsr^ they chose 1 maser spot for which they derived positions at 16 epochs during 
this 2 year period (Figure [5]l. They derive a parallax of vr = 2.29 ± 0.10 for this 
maser spot, corresponding to a distance of 437 it 19 pc. Their result does assume that 
the space velocity of the maser is constant and is not being accelerated in an outflow or 
disk. Similar to the issue for radio stars, resolvable variations in source structure could 
impact the precision of such distance determinations. 
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Figure 5. Space motion of Orion sources observed with very long baseline radio 
interferometry and used for distance determinations. Left: space motion of the non- 
thermal radio star GMR A over a 2 year period (From Sandstrom et al. 2007 1; Right: 
space motion of a water maser spot in the Kleinmann-Low Nebula over a 2.5 year 
period (From Hirota et al. 2007 1. The proper motion plus parallax reflex motion best 
fit is plotted in each case. 



2.4. Distances Based on Reddening 



Brown et al.| also considered the results of their VLU BW data in a traditional Ay vs 
distance modulus plot to estimate a distance to the Orion A Molecular cloud. In a 



finding repeated by the work of Wilson et al. ( 2005 1, Brown et al. find that the near 
edge of the cloud begins to increase the measured Ay values at a distance modulus of 



320ib70 mag. From a comparison of Ay^^ 100 //m emission from IRAS, Brown et al. 
find that the far edge of the cloud is at 500±30 pc. 

Wilson et al.| (|2005 ) report on the work of Wilson ( 2001 1 who studied the variation 



in color excess of the Hipparcos stars toward the Orion A cloud as a function of their 
Hipparcos distance. Finding the distance at which reddening increases substantially 
is interpreted as the distance to the cloud. These authors report an apparent distance 
gradient, ranging from the northern part of Orion A where the ONC is located to the 
southern filaments, which correlates with a gradient in cloud radial velocity. 



2.5. Papers Often Cited Inappropriately as Distance References 

This section addresses two issues: first the citation of distances to the ONC which have 
no traceable scientific source, and second the citation of papers in which no distance is 
quoted or in which the distance is taken directly from another source. 

A variety of distances falling within the observed scatter in measured distances 
(Table |2]) have been assumed in recent studies of the ONC. While perhaps valid as 
estimates guided by previous literature, it should be pointed out that there are no formal 
distance estimates which correspond to often quoted values such as 440 pc (DM = 8.22 
as assumed by |Herbig & Temdrup|1986 l, 450 pc (DM = 8.27, as asserted by |Genzel & 
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Figure 6. Distance estimates to Orion Id/Ic/Ic* as a function of publication date. 
Error bars are plotted where available. See Table |2] and text for description of sym- 
bols. 



Stutzki|[T989 seemingly an average of the Ic and Id distances from Warren & Hesser 



1978), 470 pc (DM = 8.36 as cited by |Genzel & Stutzki| 19891 for the result of Genzel 
et al. 1981, though not what is quoted in the original paper).These distances appear 
to be either round number estimates of bona fide published values, ad hoc averages of 
some sub-set of published values, or unpublished revisions or restatements. 

Often cited, but inappropriate, references for the distance to the ONC are the fol- 
lowing. The original paper segregating the Orion associations la,b,c,d according to 
morphology, by [Blaauw (1964), quotes a distance of 460 pc (distance modulus 8.3 
mag) which is directly referenced to a work by Borgman & Blaauw | ( |19^4| ) making the 
|Blaauwl paper not an original source. Next, the well-cited comparison of open cluster 
color-magnitude diagrams and main sequences performed by Mermilliod ( 1981|) pro- 
vides a distance to Orion with no explicit sample, methodology, or error. In general 
distances in this work are from zero age main sequence fitting; Orion is placed into a 
group with NGC 623 1 and NGC 2264 for this purpose. Finally, the extremely valuable 
proper motion study of Jones & Walker ( [1988[ ) was not, however, a distance estimate 
to the Orion Nebula Cluster. These authors simply showed that the distance of Walker 



( 1969 1 was consistent with the rejection of the number of foreground objects expected 
based on luminosity function analysis via the proper motions distribution; a similar 



analysis and conclusion was drawn by van Altena et al. ( 1988 1. 



2.6. Final Distance Thoughts 

As evidenced in Table [2j there is not only a large range in the distance estimates to the 
ONC region, but most measurements are accompanied by a large error bar, 15-20%. 
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It is interesting to illustrate the scatter in measured distances to the ONC region by 
plotting the derived distance modulus and error as a function of time, as in Figure [6] 
Those distance estimates that are the most reliable (filled symbols) are distinguished 
from the others that either are not exclusively derived from Id members or use less 
robust techniques. A notable feature of this figure is that beyond 1950 the scatter in 
distance estimates is relatively unchanged in time. While there is a clear upper bound to 
the distance measurements at ~ 525 pc, there seems to be some emerging convergence 
at 400 pc from the many varied techniques used in recent years. 

A problem that persists in more accurately separating the distance to the Orion 
Nebula Cluster from that to the Orion Ic* association is that the two subgroups are pro- 
jected along the same line of sight, where the Ic* group is concentrated primarily along 
and parallel to the Orion A Molecular cloud, which contains the slightly younger ONC 
and Trapezium (See, for example, Figure[T]l. In general, there appears to be consistency 
in the relative distances between the Orion subgroups amongst those authors quoting 
values for various of the la, lb, Ic, Id regions. Although the absolute distances have 
systematic offsets between authors and techniques, Orion Id is typically found to be the 
furthest. 

Indeed, the Id and Ic* subgroups are so aligned that a review of the membership 
statistics from the [Hille nbrand ( 1997),|Rebull|(pfel), and| Carpenter et al. |(|200]l) wide- 
field studies of this region reveals no morphological signatures that can separate the 
two entities. Distributions of infrared excess stars are more or less continuous from 
the NGC 1977 HII region down to the NGC 1999 clustering south of the ONC. Sub- 
clustering as seen by eye, and in Ha surveys are probably the result of extinction rather 
than well segregated clusters. The only significant physical difference between the Ic* 
and more embedded stars appears to be differences in their typical ages. The ages of 
these subgroups are not included in this discussion, but the historical age estimates 



among the spatially defined subgroups do seem to distinguish themselves (e.g. Blaauw 
[Warren & Hesser, ,Brown et al., ) with Id the youngest. 



3. Characterizing the T Tauri Population 
3.1. Optical to Infrared Imaging Surveys 

Photographic and Photoelectric Surveys There have been a significant number of 
photoelectric surveys published before the CCD era that include tables of source pho- 
tometry that could have utility in variability studies. Table [3] was created for the pur- 
pose of documenting these sources of literature photometry. It is supplemented with 
details from a few very large photographic surveys, e.g., [Andrews & GrossTe ([1981 1 



with Ni, ~ 15000. Data from some papers, such as that data listed in Warren & Hesser 



(1977)'s primary UBV table, are a merger of a very large number of published data 
sources (in that case 35 separate papers), and the references in these amalgamations are 
generally not reproduced here. 

Ha Surveys Slitless optical grism surveys of young star forming regions can be valu- 
able tools for identifying young stars (e.g. Herbig & Bell|1988 l. This is because young 



stars frequently show strong Ha line emission, which is related to their active chro- 
mospheres as well as circumstellar accretion; however, the very strong hydrogen line 
emission background of the Orion Nebula probably result in a significant underestimate 



of the true membership if based upon Ha statistics alone. Haro ( 1953j ) documented 255 
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Table 3. Large scale photographic (pg) and photoelectric (pe) data sources for 
the Orion Nebula. Ordered chronologically by epoch of publication. No epoch is 
recorded when no detailed epoch of observation could be determined or the data 
were the average of many other sources. 



Paper Data Epoch Filters A^* 



Parenago'('1954a| 




pg 




rUpg 


2982 


bharpless ( ly^^l 




pe 


1951 


BV 


190 


ISharpless ( iy5?| 




pe 


1951 


UB 


184 


Johnson ( l957r 




pe 


1954-1955 


UBV 


49 


Sharpless ( iy62| 




pe 




UBV 


~ 180 


Kopylov & Stralzys'('1963') 


pg 


1952? 


nipg 


- 2000 


Morgan & Loden (1966 J 


pe 




UBV 


36 


Lee (196«| 




pe 




UBVRIJKL 


196 


Walker ( 1969|l 




pe 


1958-1967 


UBV 


- 300 


Penston (19 /3 1 




pe 


1970-1971 


UBVRIJHKL 


51 


Penston efaTT 19751) 


pe 


1972-1973 


UBVRIHKL 


48 


McNamai-a (1976 


bF 


pe 




UBVRIJHKL 


51 


Warren & Hesser 


(jl9/ 7 1 


pe 


1968 


UBV 


526 






pe 


UBV 


109 






pe 


1972 


uvby 


492 


Shevchenko'('1981| 


pe 


1978-1979 


UBVRI 


117 


Andrews & Grossiel(|1981|l 


Pg 


1979 


UBVI 


~ 15000 


Isobe & Sesaki (1982) 


pg 


1970-1971 


RI 


413 


Rydgren& Vrba,(,1984|) 


pe 


1982-1983 


UBVRI 


41 



Ha stars within a 3.5 degree region around the Trapezium while IParsamian & Chavira 



( 1982) cataloged 534 Ha stars in a 5 degree region. The Kiso Orion surveys (Wirami- 
hardja et al.|[T989t |Kogure et al.1[T989, ,Wiramihardja et al.|[T99T| [19931 [Nakano et al. 
1995| ), while valuable for their coverage of most of the Orion constellation, appear to 
be very incomplete in the ONC as evidenced by the lack of a strong peak in the stel- 



lar density within the Nebula as found by the subsequent analysis of Gomez & Lada 
( 1998 ); Jones & Walker (see also 1988| ). In their review of the Orion association, Brand 
& Wouterloot ( 1991 ) collated the existing Ha star catalogs into a single list, including 



87 new stars from Wouterloot & Brand ( 1992 ). 



Modern Optical Surveys A review of modem optical CCD surveys of the Orion Neb- 



ula begins with the work of Herbig & Terndrup ( 1986 ). Their CCD observations were 
taken with the 40 inch Nickel telescope at Lick Observatory, had a pixel resolution of 
0.267" and consisted of a mosaic of small (~ 2.5') fields. The authors used narrow- 
band interference filters to minimize nebular contamination but final photometry was 
reduced to and reported in the Johnson-Cousins Vic system. Their Table 1 contains 



photometry for 98 of the 140 sources detected and uses the Parenago (1954a) number 
system except for 30 sources that are labeled "anonymous." The authors used these 
new data to construct color-magnitude diagrams and explore the age and age spread for 
the cluster, finding most stars to be ~ 1 Myr or younger (Section 3.4. | ). 

The first ONC photometry from the Hubble Space Telescope was published by 



Prosser et al. ( 1994 ). This survey consisted of 1 1 irregularly mosaicked Planetary Cam- 
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Figure 7. The Orion Nebula in optical light as mapped by the Advanced Camera 
for Surveys with the Hubble Space Telescope. The field is approximately 30' x 30' 
in size with North oriented up and East to the left. From the Orion HST Treasury 
program, PL M. Robberto. 



era fields in the F547M and F875M filters. Their Table 4 contains aperture photometry 
for 326 objects, using an aperture beam of 0.12" and converted into the Vic passband 
system. Cross references from their "PC" identifier system to that of Jones & Walker 
(JW) and Parenago (P) are given. Unfortunately, they report their astrometry to be 
quite poor {r^ 1"). The high resolution of these data provided excellent new statis- 
tics on visual binaries in the cluster, identifying 35 sub-arcsecond pairs (their Table 6). 



Additional HST observations were obtained, reduced and analyzed by Robberto et al. 



(2004 1. Their results include observations in the F336W, F439W filter passbands as 
well as data from archived F547M, F791W images. They tabulate the resulting U BVI 
data for 40 sources with spectral types from |Hillenbrand ( |1997 1. 
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The comprehensive survey of the Orion Nebula Cluster by Hillenbrand (19971 
included new Vic photometry in addition to a large corpus of spectroscopy (Section 
3.3. 1. Hillenbrand cataloged 1578 sources including 332 new detections with approx- 
imate completeness limits of Ic ~ 17.5 and V ~ 19. Her tabulation was constructed 
from new data (3 epochs) and literature sources (7); photometry of sources appearing 
in multiple catalogs was chosen based on the angular resolution of the survey (e.g., the 
Prosser et al. HST results were given preference). Their numbering system is a merger 
of Jones & Walker (JW; #1-1053), Parenago (32 sources, e.g. Parenago 1891 = Hillen- 
brand 1891 = 6*1 Ori C), Prosser et al. HST sources (9000-i-PC#) and new detections: 
3000-i-N for epoch 1993 data, 5000-i-N for epoch 1995 data and 6000H-N for epoch 1996 
data. Most of the global stellar properties of the Orion Nebula members are derived 
from this comprehensive study. 

In addition to these published surveys modem telescope archives contain large 
quantities of publicly available optical data. The most significant of these is the 104 or- 
bit Cycle 13 Hubble Space Telescope Treasury Program (PID 10246; PI. M. Robberto) 
that surveyed a ~ 20' x 20' region of the Orion Nebula with the Advanced Camera for 
Surveys (ACS). Observations took place between October 2004 and April 2005 and the 
surveyed filters included F435W, F555W, F658N, F775W and F850LP Parallel obser- 
vations were also obtained with the Wide-field Planetary Camera 2 (F336W, F439W, 
F656N, F814W) and NICMOS (NIC3: Fl lOW, F160W); all these data can be obtained 
via the Multimission Archive at STScI. Figure [7] is from their ACS mosaic (Press re- 
lease STScI-PR06-01). A similarly large set of ground-based optical data (PID 273.C- 
5042(A)) is publicly available from the ESO archive. It was observed during January 
2005 with UBVIc and Ha filters (Da Rio et al., in preparation). 



Optical Variability Surveys Occasionally, multi-epoch variability surveys publish cal- 
ibrated time-averaged photometry for their sources. Sources for such photometry in 
the Orion Nebula include Stassun et al.| \\999\ and [Herbst et aL] ( |2002| ). The former 
provides data for their 254 periodic stars in the Ic passband while the latter provides 
narrowband photometry for 1562 objects time averaged over 45 days. The narrowband 
filter used by Herbst et al. was centered at 815.9 nm. In both cases the reported peak to 
peak variation of 0.2 magnitudes is probably a good measure of the typical uncertainty 
inherent in a single epoch optical survey of a young cluster. Note that [Herbst et al. 
p002) adopted the same numbering system of Hillenbrand ( 1997| l in their Table 1, ex- 
tending it to N> 10000. The optical photometry for variable and periodic stars surveyed 
in the outer ONC by 'Rebull et al. (2000) and Rebull|(|2001|) is single epoch. 



,Rebull etal.lpOOO^ and ,Rebull,(2001^ presented UVI photometry of the "flanking 
fields" of the ONC. The RebuU flanking field surveys cover an area from 1.5 > i? > 
0.25 degree from 6^ Ori C out to the limits of the Orion A cloud to the north (see chapter 
by Peterson & Megeath) and south of l Ori and the OMC-4 clumps. The Tables 1 and 
2 of [RebuU et al.| ( |2000| ) include a total of 4792 (1620) sources (candidate members) 
that were found in or near the locus of confirmed Orion sources on the optical color- 
magnitude diagram. A total of 1564 sources (726 candidate members) have valuable 



wide-field U band data in Rebull et al. (20001. 



Near-Infrared Data Near-Infrared observations are necessary to explore embedded 
populations in young regions like the Orion Nebula and are sensitive to re-radiated ther- 



mal emission from circumstellar disks. Penston ( 1973 1 performed an optical-i-infrared 



survey of 5 1 Parenago stars over a 0.5° region. These observations included the first and 
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essentially only wide-field 3 /^m photometry for Orion Nebula members until the year 



2000. Additional HKL data for 35 stars were presented in Pension et al. ( 1975 1. Lons- 



dale et aiy] ( |1982[ l expanded the census of sources near the B-N object. Although they 
tabulate only sources within 35" of the B-N object, larger maps including the Trapez- 
ium are referenced and shown, including a source map from Becklin et al. ( 1976[), who 



surveyed at 2 and 20 /im but did not tabulate any point source photometry. Hyland 



et al. ( 1984*) produced a non-chopped 2 /im map that covered the entire OMC-1 molec 



ular ridge (both B-N/K-L and OMC-IS) and tabulated 88 sources with cross-references 
to Lonsdale et al. and Parenago. 

A cruc ial near-infrared survey of the Orion Nebula was performed by [McCaugb 



rean & Stauffer (1994). These authors obtained two complementary sets of K' data, 



covering a total of 82'' x 82" and centered on the Trapezium stars. The higher reso- 
lution tip-tilt corrected images had a final spatial resolution of 0.35" and the authors 
quote an astrometric precision of 0.06". Their Table 1 lists 123 detections (48 new 
stars), including photometry even for the brightest OB stars and extending to a quoted 
completeness limit of K' = 16. This tabulation is the origin of the "TCC" or Trapez- 
ium Cluster Catalog identifiers and provides cross-references of their near-IR data to the 
catalogs of Jones & Walker (1988), Parenago (1954a) and Prosser et al. (1994) sources, 
as well as a detailed and valuable cross-referencing of known VLA radio sources from 
Felh et al.|(|1993b| ), the proplyds (e.g., |0'DeU & Wen.1994^, Ha sour ces from ILaques] 



& Vidal (19791, mid-IR sources and structures from Hayward et al.| ([r994) and their 



Table 1 of optical/near-IR sources. 

Subsequent near-IR surveys can be divided into those which covered very large 
areas of the Nebula and those that concentrated on the central ~ 5' around the Trapez- 
ium. Wider field surveys that provide large tabulations of near-IR photometry include 



Ali & Depoy ( 1995 1, Hillenbrand et al. ( 1998 1 (multi-epoch) and Carpenter et al. ( 2001 1 



(time series). Concerns about questionable and probably nebular extractions should be 
applied to the results of almost any survey for sources against the bright, background 
Nebula. [Hillenbrand et al.| ( |1998| ) and [Hillenbrand & Carpenter] ( |2000j ) suggested a 
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large percentage of false and duplicate detections in the Ali & Depoy data, while Hil- 
lenbrand & Carpenter also noted 1 or 2 nebular knots among the sources listed in the 
optical catalog of Hillenbrand (1997). Similarly, in a 17 x 17' field centered on the 
Trapezium, Getman et al. ( 2005c| ) found 1145 sources in the 2MASS point-source cat- 
alog that lacked Chandra X-ray detections; however, only ~ 200 of these are good 
quality 2MASS detections; the rest (~ 900) lack detection in multiple bands and most 
are probably spurious. The catalog of [Hillenbrand et aL ( 1998[ ) avoids this problem 
because it tabulates near-IR photometry only for those 1578 optical sources listed in 
[Hillenbrand, ( , 1 997 ) ; this tabulation includes new bolometer and array observations sup- 
plemented by literature results. Similarly, the Carpenter et al. (2001 1 2MASS near- 
infrared variability survey tabulates photometric results for those 1235 variables (out 
of 17,808 sources) found in a 0.84° x 6° region. Variable stars in this catalog were 
typically observed 16 times over a 2 year period. 

Narrow field, deeper surveys of the region immediately around the Trapezium 
have i ncluded the K-only AO survey of Simon et al. (1999), the Keck ( HK) sm- 



vey by [Hillenbrand & Carpenter] ( |2000| ), th e HST-NICMQSX/-g) survey by [Luhman 
etal] ( [2000| l, UKIRT IJH observat i ons by [Lucas & Roche| ( [2000] ), a mu lti-telescope 



JHKKsL survey by Muench et al.| (2002!) and a Gemini JHK survey by Lucas et al. 
(2005 1. A comparison of the regions surveyed by most of these authors was given by 



Muench et al. ( 2002[ ). Additional near-IR data are included in the COUP catalogs and 



include previously unpublished photometry derived by McCaughrean. As mentioned 
previously, archival NICMOS observations obtained in parallel to ACS imaging of the 
Nebula provide a non-contiguous but as yet unpublished data set for future use. Addi- 
tionally, the CFHT archive contains a large WIRCAM/UKIRT YJHKsH2 data set of 
a large field surrounding the Orion Nebula. 



Thermal and Mid-Infrared Data Initial mid-IR scale maps were limited by a combi- 
nation of the very strong nebular background and the poor angular resolution of early 
mid-IR cameras. Works by [Ney & Allen[ ( [1969[ ); [Lemke et al.[ ( [1974[ ); [Fazio et al 
( 1974); Gehrz et al. (1975]), which span the wavelength regime from 20 to 100 ^m, tell 
us little about the overall stellar content or properties of the members of the star cluster 
although they do reveal many details about the structure of the photodissociation re- 
gion, e.g. the modem study by Kassis et al."(2006). Even today the 20 ^m flux from the 
Nebula overwhelms all but a few bright protostars and, for example, renders the inner 
15' of the Nebula saturated at 24 fiva with the Spitzer Space Telescope. 



There is a multitude of mid-IR surveys (e.g., Rieke et al. 



1978 



2004 



Lee et al.[[T983, ,WynrrWimams et al.|1984, ,Gezari et a 



Rieke et al. 


1973; Beichman et al. 


, ,Gezari et a 




1998 


Shuping et al. 



Greenhill et al. 2004a i that have focused only on small embedded regions Uke the 
B-N/K-L. Full discussion of the stellar content of these embedded regions is reserved 
for a subsequent chapter (O'Dell et al.; Part II). The 12 //m image of the Ney- Allen 
region ( Ney et al.|1973| ) of the central Trapezium from McCaughrean & Gezari ( 19901 
was the first modern mid-IR array observation, revealing many narrow arcs and details 



of the prominent structures in this region. Hay ward et al.[ ( [1994, ) published 8.8 and 
11.7 //m SpectroCam (SC) images of the central region, including fluxes for 13 sources. 
Hayward ( |1994 ) present additional 10 and 20 fim maps of the central Trapezium. 

Thermal infrared 3 fim data for ~ 400 sources in a 7' x 7' region were analyzed 



by |Lada et al. ( 2000| ), who used these data to estimate the disk fraction as a function of 
source mass and to identify a large sample of protostars throughout the OMC-1 cloud. 
The data tables used in that work were presented in [Muench et al.| ([2002|). Deeper, 
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Figure 9. A comparison of the Trapezium at mid-IR (11.7 /zm, left) and optical 
(HST WFPC2) wavelengths. Figure reproduced from |Smith et al.| ( |2005| l. It is im- 
portant to note that at wavelengths longer than these, e.g., > 20 /im, the only promi- 
nent point sources that can be detected against the strong nebula are located near the 
B-N/K-L region (|Robberto et al.|2005[). 



higher resolution 3 fim data of ~ 400 sources in a smaller 5' x 4' region were published 



by Lada et al. ( 2004^ ); these data extended measurements of the disk fraction into the 
brown dwarf regime and provided new results on the protostellar population. 

Recent longer wavelength observations with the spatial resolution necessary to 



detect individual sources against the bright background include Rob berto et al. (2005 1 
(10 & 20 /um; 4' x 5'; 0.5" resolution; 177 sources) and [Smith etal . (2005) (11.7^m; 
2' X 3'; 0.35" resolution; 91 sources), see Figure [9] These works focused their study 
on the proplyds, jets and emission structures in the PDR, which are a focus of the 
following chapter (O'Dell et al., Part II). Publications that include Spitzer observations, 
which are much lower resolution than any of the ground based mid-IR observations, 
include 'Rebull et aL| (|2006[) and|Cieza & Baliberl(|2007'); both of these works focused 



on the relationship between the rotational properties of the young stars and their disk 



excess properties (Sect. 3.6 



3.2. X-ray Observations of the Nebula 

The birth of stars takes place in thermodynamically cold and neutral media with charac- 
teristic energies of <C 1 eV per particle. Paradoxically, those processes associated with 
star formation produce and are subject to violent high energy processes with character- 
istic energies of ^ 10'^ eV. The principal evidence for this is X-ray emission from stars 
throughout their pre-main sequence (PMS) evolution. The ONC was the first cluster 
of PMS stars to be detected in the X-ray band ( Giacconi et al.|1972 i and non- imaging 
studies soon found that the X-ray emission is extended on scales of a parsec or larger 
(Bradt & Kelleyf 1979) 1. Early explanations for the Orion X-rays included winds from 
the massive Trapezium stars colliding with each other or the molecular cloud, and hot 



corona or magnetic activity in lower mass T Tauri stars. The Einstein ( Ku & Chanan 
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Figure 10. The centr al 4' x 4' COUP A CIS-I field shown with 0.25" pixels. Im- 
age reproduced from Get man e t al.|( 2005c| . The displayed image has been smoothed 
from the natural integer version using the Chandra Interactive Analysis of Observa- 
tion (ClAO) csmooth procedure. Red, green and blue colors correspond to counts 
in the 0.5 1.7 keV, the 1.7 2.8 keV, and the 2.8 8.0 keV bands, respectively. 



19791 ), ROSAT (Gagn e et al. 11995] ) and ASCA ( |Yamauchi et al.|1996| ) imaging X-ray 



observatories established that both the massive Trapezium stars and many lower-mass 
T Tauri stars contribute to the X-ray emission. 

The ONC was intensively studied during the first year of the Chandra mission 
with several instrumental setups: the Advanced CCD Imaging Spectrometer (ACIS) in 
imaging mode ( [Garmire et al. 2000 Feigelson et al. ^002a b , 2003]) and as detector 



for the High Energy Transmission Grating Spectro meter (|Schulz et al. 2000 2001 
2003a|b[ ), and with the High Resolution Imager (HRI; [Flaccomio et al.|2003a|b] ). While 



many valuable results emerged from these early Chandra studies, it was recognized 
that more would be learned from a deeper and longer exposure of the Orion Nebula 
region. During the fourth year of its mission Chandra performed an unprecedented 
10 day (net exposure) nearly-continuous observation of the Orion Nebula, nicknamed 
the Chandra Orion Ultradeep Project (COUP). 



Chandra Orion Ultradeep Project The COUP study detected more than 1600 X-ray 
sources, ~ 1400 of which are young stellar objects (Get man et al.||2003c ). Figure 10 
shows a "true-color" X-ray image of the central 4' x 4' Trapezium region centered on 
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the larger 17' x 17' field of Chandra-AClS-l. Absorbed COUP sources appear blue 
and unabsorbed sources appear red. Chandra X-ray studies are particularly effective 
in uncovering heavily obscured low-mass cloud populations (X-rays penetrate up to 
hundreds of magnitude of absorption into the cloud) and in discriminating cloud PMS 
populations from unrelated older stars (X-ray emission from PMS stars is 10^ — 10^ 
times elevated above main sequence (MS) levels). Most of the non-PMS contaminants 
in the COUP field are extragalactic active galactic nuclei (AGNs), which can be con- 
fused with non-flaring YSOs. (Only 16 probable field stars with discrepant proper mo- 
tions and NIR colors are present in the COUP source list, which are available through 
Getman et al.||2003b]). But the long exposure improves the opportunity for captur- 



ing powerful X-ray flares which are characteristic of YSOs and not AGN. Based on 
the variability analysis of heavily absorbed COUP sources without optical/NIR coun- 



terparts and detailed simulations of the extragalactic background population, Getman 



et al. (2005a I argue that 75 COUP sources are previously unknown embedded cloud 
members, of which forty-two are confirmed by the detection of powerful X-ray flares. 
These X-ray discovered stars are spatially clustered within the two well-known OMC- 
1 cores and the dense molecular filament, which extends northwards from OMC-1 to 
OMC-2/3. 

The census of COUP sources with confirmed Orion membership includes 1315 
stars with known optical/NIR counterparts, 75 new embedded stars, 16 unidentified 
likely new lightly obscured members of ONC ( Getman et al.|2005c i, and two faint X- 
ray sources associated with the Herbig-Haro object HH-210 (COUP #703 and #704; 
Grosso et al.||2006 l. Three classes contribute roughly equally to the integrated X-ray 
luminosity in the hard 2 — 8 keV band: 10 unobscured hot Trapezium stars earlier 
than B4, 839 cool (later than B4) hghtly-obscured COUP sour ces with log Nh <^ 22 



cm (Ay ^5 — 6 mag), and 559 heavily-obscured stars (Feigelson et al. 20051. 
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Figure 1 1 . Diagram of the Orion Nebula field showing 1408 COUP X-ray sources 
assoc iated with the Orion N ebula segregated by obscuration. Figures reproduced 



fromJFeigelson et al. (2005V Left: Lightly obscured subsample with <^ 10^ 



Right: Heavily absorbed subsample with >^ 10^^ cm^^. The large 
triangles show 10 hot 07-B3 stars, while dots show the remaining cool member 
population. 
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Figure 12. Histogram with the differential distribution of hard (2.0 — 8.0) keV 
band X-ray luminosity corrected for absorption. Figure reproduced from Feigelson 
et al. (2005). Solid line denotes the COUP unobscured cool star sample with Gaus- 



sian fit, dashed line show the obscured sample, and the black-filled histogram show 
the 10 hot stars, including 9^ Ori C. 



The spatial source distribution for the cool unobscured (Figure 11 left), and heavily 
obscured populations (Figure 1 1 right) show a spatial asymmetry - a deficit of stars to 
the east on 0.5 — 1 pc scales - consistent with violent relaxation in the stellar dynamics 
(see, however, Furesz et al.||2008| ) and the concentration of obscured sources around 
both OMC-1 molecular cores. The X-ray luminosity function (XLF) of the unobscured 
cool population is > 90% complete down to M ~ O.IMq and ~ 50% complete down 
to M ~ O.OSMq. The XLF shape is roughly log-normal in shape and the obscured 
population is deficient in lower-luminosity stars due to localized circumstellar material 
(Figure 12 1. One-third of the Orion Nebula region hard-band emission is produced by 
the bright O star 6^ Ori C, and half is produced by lower mass pre-main sequence stars 
with masses 0.3 < M < 3 Mq ( [Feigelson et al.|2005| ) 



With the detection limit of Lx,min ~ 10^*^ ergs/s for the unobscured COUP pop- 
ulation. X-ray emission was detected from more than 97% of the optically visible late- 
type (spectral types F-M) T Tauri stars (TTS) in the ONC, demonstrating that there is no 



'X-ray quiet" population of late-type stars with suppressed magnetic activity. Preibisch 



et al. ( 2005a[ ) show that TTS with known rotation periods lie in the saturated or super- 



saturated regime in a diagram comparing X-ray activity with the stellar interior Rossby 
number (Figure 13 left). But the TTS show much larger scatter in X-ray activity than 
main sequence stars. This scatter is partly attributable to accretion: while the X-ray ac- 
tivity of the non-accreting TTS is consistent with that of rapidly rotating MS stars, the 
accreting stars are less X-ray active (by factors of ~ 2 — 3), perhaps because magnetic 
reconnection cannot heat the dense plasma in mass-loaded accreting field lines to X-ray 
temperatures. The fact that COUP stars do not show the drop-off in magnetic activity as 
stars rotate more slowly may suggest that the magnetic dynamo process is saturated in 
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Figure 1 3 . Left: Fractional X-ray luminosity versus Rossby number for the COUP 
T Tauri stars. Figure reproduced from Preibisch et al. ( 2005a[ l. The grey shaded area 
shows the relation and the width of its typical scatter found for MS stars; Right: The 
black line shows the integrated distribution of X-ray luminosities for COUP Orion 
stars (ages 0.1 — 10 Myr) in a narrow mass range. Comparison distributions are 
the Pleiades cluster (red line, age ^ 100 Myr), Hyades cluster (blue line, age ^ 500 
Myr), and solar neighborhood stars (green line, ages ~ 1 — 5 Gy). It is clearly shown 
that the activity-age relation continues through the pre-main sequence phases. Figure 
reproduced from Preibisch & Feigelson (2005|l. 



some way and/or that a different dynamo is operative in young stars that is independent 



of rotation. Preibisch et al. (2005a ) do find that COUP X-ray luminosities are correlated 



with stellar mass and volume, which generally suggests a turbulent convective dynamo 
model. 

For main sequence stars older than 50 Myr, it has long been known that younger 
stars are more magnetically active than older stars. Preibisch & Feigelson ( ,2005) clearly 
establishes that the activity-age relation continues through the PMS phases (Figure 



13 right) and find a decay law that is mass-dependent at young ages. Wolk et al. (2005 1 
used a complete sample of 1 solar mass Orion stars in the COUP field to show that 
analogs of the young Sun spend one-fourth of their time in flare state, exhibit incredibly 
high levels of magnetic activity with the median luminosities 2-3 orders of magnitude 
higher for both "quiescent" and peak flare levels compared to the contemporary Sun. 



Caramazza et al. ( 2007 1 further show that X-ray flare frequency in young lower- mass 
(0.1 — 0.3Mq) stars is indistinguishable from that of the young solar analogs. Finally, 
Maggio et al. ( 2007] ) find that elemental coronal abundances in X-ray luminous young 
Orion stars are similar to those of older magnetically active stars. 



X-rays from Embedded Sources The spatial distribution of "obscured" COUP sources 
clearly traces the basic structures of the central cluster; each of Trapezium core, the 
B-N/K-L and OMC-IS regions appear as over-densities in Figure [T4| The detailed 



properties of the COUP detected X-ray embedded sources in B-N/K-L and OMC-IS 
regions (see boxes in Figure 14 1 are discussed by |Grosso et al. (20051. Grosso et al. 
found 60 COUP X-ray sources toward the OMC-IS dust continuum core, with more 
than 60% of them being obscured. In the B-N/K-L region 43 sources were detected and 
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Figure 14. Distribution of COUP sources differentiated by X-ray absorption col- 



umn density and compared to the SCUBA 450 /im map of [Johnstone & Bally| (|1999 
shaded contours). Black dots: obscured (A'h > 10^^ c 



cm~^) sources; white dots: 
Three major sub-clusterings 
are differentiated: ordered from North to South these include t he B-N/K-L (boxed) , 
the Trapezium, and OMC-IS (boxed). Figure reproduced from Grosso et al. (2005 1. 



lightly-absorbed (< 10^^ cm ^) COUP X-ray sources 



half of these were obscured. Based on comparison of the X-ray luminosity function of 
the observed X-ray populations embedded in OMC-IS and B-N/K-L with that of the 
unobscured ONC population, , Grosso et al. estimated total populations of 70 versus 80 
embedded stars residing inside the OMC-IS and B-N/K-L cloud cores, respectively. 

Close-up images of these two regions, scaled to the same physical size are com- 
pared in Figure 15 This is the first direct measurement of the low-mass population of 
the OMC-IS cluster with 18 new X-ray sources without infrared counterparts. COUP 
OMCl-S detections include the most embedded X-ray source in the COUP survey, 
COUP 632 (= TPSC 1), a protostar with Ay ~ 500 mag of visual absorption. X-ray 
sources are found close to four luminous mid-IR sources - B-N, IRc3-i2, IRc2-C, and 
Source n - but their X-ray variability and spectral properties are typical of coronal ac- 
tivity of low-mass companions rather than wind emission from massive stars. No X-ray 
emission is seen from the radio-bright massive protostar Source I. 

Using the combination of compiled CTIO-ISPI near-IR with Spitzer IRAC mid- 
IR imaging data, Prisinzano et al. ( |2007 1 establish the list of 45 protostellar candidates 
within the COUP field of view: 23 designated as Class 0-Ia with their IR SEDs mono- 
tonically rising from to 8 fim and 22 designated as Class O-I with SEDs rising from 
K up to 4.5 /Lim. Of these, 62% have X-ray counterparts in COUP data. Their tabu- 
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Figure 15. A comparison of two embedded subclusterings in the Nebula as traced 
by COUP multi-spectral data. The figures show the B-N/K-L (left) and OMC-IS 
(right) subclusterings, corresponding to the regions outlined in Figure 14 Color 



coding (red, green, and blue) correspond to photons in the 0.5-1.7 keV, 1.7-2.8 keV, 
and 2.8-8.0 keV bands; red markings correspon d to well known prot ostars, hot cores 
or far-infrared source. Figures reproduced from Grosso et al. (2005 1. 



lations also contain cross-references to a number of thermal IR surveys of the nebula. 
The spatial distribution of these protostellar candidates trace the dense molecular fil- 
ament that extends northward from OMC-1 to OMC-2/3 clouds and is similar to that 



of 75 likely new embedded cloud members found in Getman et al. (2005c I. However, 
due to nebular contamination and crowding in mid-IR, Prisinzano et al. were not able 
to classify many X-ray embedded sources located in B-N/K-L and OMC-IS regions 
( Grosso et al.|20 05). In addition, a sub-cluster of seven highly embedded X-ray sources 
in OMC-IS (COUP # 582, 594, 615, 633, 641, 659 and 667) at approximately RA,DEC 



= 05:35:14.8, -05:24:12 (J2000; see also Figure 15 1 is simply invisible in Spitzer data. 
In regards to the evolution of the X-ray emission, Prisinzano et al. find that Class 0-Ia 
protostellar candidates are intrinsically less luminous than the Class II stars. 



X-ray Properties of Orion Brown Dwarfs In the core of the ONC, 9 spectroscopically- 
identified brown dwarfs were detected with the COUP ( jPreibisch et aL|2005 bi). The low 
detection rate is seemingly related to the substantial extinction to most of these brown 
dwarfs (Figure [T6]). There is no evidence for changes in the magnetic activity around the 
stellar/substellar boundary; instead the X-ray properties of the detected brown dwarfs, 
including spectra, fractional luminosities, and flare rates, are similar to those of the low- 
mass stars in the ONC. Trends in the fractional X-ray luminosity and surface flux and a 



comparison to the X-ray properties of late type field dwarfs led Preibisch et al. ( 2005b I 
to conclude that the photospheric temperature of late type stars rather than source mass 
or surface gravity controls the X-ray emission mechanism. 



X-ray Observations of Orion Flanking Fields Because the COUP survey samples 
only the inner parts of the Orion Nebula, X-ray surveys of what have been deemed 
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Figure 16. The relationship between optical extinction and X-ray luminosity (or 
its upper limit as marked by arrows) for 34 Orion brown dwarfs. Figure reproduced 
from Preibisch et al. (2005b|. Detectability of Orion brown dwarfs in X-rays by the 
COUP was limited, apparently, by extinction. 



the Orion "Flanking fields" are important. This is because they provide membership 



information at large cluster radii, where other methods are ambiguous. The Ramirez 



et al. ( 2004 1 survey follows the axis of the cloud, sampling Chandra ACIS fields north 
and south of the immediate Nebula. While they are much less sensitive than the COUP 
observations, their results provide interesting evidence for age gradients along the axis 
of the cloud. New surveys with XMM covering all of the OMC and Orion A clouds are 
currently being completed and will soon yield more spatially complete results as well 
as extensive publicly available archival data. 



3.3. Spectroscopic Surveys 

Historic Studies The "Orion population" stars have long been of interest to spectro- 
scopists. Early and numerous studies using objective prism plates were published by 
e.g. Herbig, Haro, Sharpless, Johnson, Walker, and Pension among others. These au- 
thors were interested in emission-line aspects as well as spectral types over the Orion 
Ic and Id regions. Specifically, spectral types for stars in the ONC, as defined above, 
were reported by: Blanco ( 1963 ); [Co hen & Kuhi|(|1979|l;|Greenstein & Stru ve ( 1946^ 
Herbig as quoted infmIEr|(1969); Johnson (1965); Levato & Abt ( 1976); Abt & Le 



vato 



([19771); [Lallemand et al, ( ,1960^; iMcNamara, ( |1976a| ); [Parenago, ( , 1954a| lT [Penston 



([1973 ); Pension et al. ( 1975 1; Strand ( 1958 ) though referenced mostly to Sharpless; 



Smith et al.[ ( |1983| ); |Trumpler[ ( [1931| ); [Walker[ ( [1983| ). Many of these studies also in- 
clude the larger Orion population. In the majority of cases the literature of this era was 
focused on identifying the emission line stars, on characterizing the cluster sample, 
and on understanding whether the fainter objects should be interpreted as reddened, or 
lower mass, or of older age. 
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Modern Low Resolution Surveys Modem optical spectra of a few tens of ONC stars 
were produced by pioneers such as Herbig & Temdrup] ( 1986) and van Altena et al 



([T988) , while photographic Ha plate surveys have also continued (see previous Sect. 
3.1.[). The largest collection of published spectral types is contained in Hillenbrand 



( 1997 1, which incorporated new data on many hundreds of stars as well as previ- 



ously published (those references above plus more modern additions from e.g. Duncan 
( [T9951 ); [Edwards et al.| ([T993]); [Wolff eTal] ( |2004l ), and unpubhshed (e.g. Prosser & 
Stauffer spectra; Samuel ( 1993[ ), PhD thesis, and Hamilton 1994, MSc thesis) informa- 
tion. Approximately 950 spectral types were provided. Since the ^Hillenbrand, (1997] ) 
publication, however, approximately 800 more spectral types over the same projected 
area have become available; an updated catalog is being prepared (Hillenbrand et al.). 

Relevant sources of new optical spectroscopy include ,Riddick et al. (2007) and 
spectral types for a few of the sources in Slesnick et al. (2004). Infrared spectroscopy 
includes Luhman et al. (20001, but has more recently focussed on the lowest mass 
candidate members of the ONC with contributions by [Lucas et"aL] ([200 1[ [2006|), and 



Slesnick et~aLl ( [2004[ ). There is also ongoing work of Lada with the FLAMINGOS 
multi-object spectrograph. Very late M and perhaps even LO or LI objects have now 
been identified in the ONC region. 



Modern Echelle Surveys Ushering in modern high dispersion studies of the Orion 
Ic/Id region were the works of [Smith et al.[ ( [TM3]l,[Walker[([l990[ ), [McNamara[ ( |T990l l, 
Abt et al.[([1991), and Duncan et al.' (1991); 'Duncan' (1993) which all focused on ro- 



tational velocitie s. King, (1993; , .Duncan & RebuU, (1996) . .Palla et al.. (2005 ) (and 
Palla et aL][2007[ ) subsequently studied lithium abundances as well as radial velocities 



for small samples of Orion stars; Sicilia-Aguilar et al. (20051 derived the same for a 
larger sample of several hundred stars, and also provided rotational velocities. Addi- 
tional rotational velocities come from Rhode et al. ( 200 1|| and Wolff et al. ( 2004^ who 



both studied stellar angular momentum. Padgett (1996 1, Cunha et al. (19951 1998), and 



Cunha & Smith (2005) published work on abundances including several ONC stars. 



Furesz et al. ( 2008 1 provide radial velocities for a large sample of ONC stars in a study 



of cluster kinematics. Most recently, Yang & Johns-KruU ( 2007 1 have achieved the 
means to study the magnetic fields of ONC stars. 

At this time, there is a considerable amount of data across the stellar mass spectrum 
on rotational velocities in the ONC that, together with rotation periods, are providing in- 
sights into stellar angular momentum at the earliest stages of stellar evolution. Lithium 
samples are far smaller, but remain valuable. Abundance information and magnetic 
field measurements are intriguing but remain rather limited at present. 



3.4. Stellar Properties 

The photometric and spectroscopic surveys described in the previous two sections have 
provided the panchromatic data necessary to assess the stellar population of the ONC. 
Compared to other such studies of older clusters, significant challenges are posed in 
the ONC by the effects on observed colors of nebular contamination and disk accre- 
tion (both of which cause blueing), and circumstellar and interstellar dust (reddening). 
Spectral continuum and certain spectral lines are also affected. Nevertheless, perse- 
verance has led to understanding of typical ages and to estimates for individual stellar 
masses. 
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Figure 17. Low mass Hertzsprung-Russell Diagram for the Orion Nebula Cluster 
(this is an updated version of that constructed in Hillenbrand' 1997). Isochrones and 
fixed mass evolution tracks are from |D'Antona & Mazzitelli, ( ,1997) . 



Stellar Ages and Star Forming History Notwithstanding subsequent criticism of his 
data quality by Walker (1956), Parenago ( 1954a I was first to observe a "cloud" of sub- 
giants later than spectral type AO and more luminous than the main sequence in the 
vicinity of the Orion Nebula. Along a pathway similar to his previous studies of NGC 
2264, NGC 6530, NGC 6611 and IC 5 146, .Walker. ( ,1969; used photoelectric photom- 
etry of sources listed in Brunf s catalog to confirm the existence of pre-main sequence 
stars in the Nebula. Other photometric + spectroscopic surveys of the 1970's followed 
suit. It was Herbig & Terndrup. ( ,1986; who quantitatively demonstrated the youth of the 
ONC stars by comparing their dereddened color-magnitude diagram to the theoretical 
evolutionary tracks of Vandenberg & Bell ( 1985| ). They found that the vast majority of 
stars were younger than 1 Myr, and therefore seemingly inconsistent in their age dis- 
tribution with the canonical 10 Myr life times of molecular clouds ( [Blitz & Shu|1980 l. 
Improvements in the data and in the theory resulted in similar conclusions being drawn 
by subsequent workers, e.g. Prosser et al.|(|1994[),| Hillenbrand| ( [T997| ) and |Slesnick et al. 
([2004). Figure [TTJis an update of the H illenbrandFi 1997 1 HR diagram illustrating this 
star forming history. Hillenbrand ( 1997| ) further suggested a mild age gradient, featur- 
ing young median ages closer to the cluster core and a slightly older population in the 
outer nebula region about 2.5 pc from the core. 

However, these later studies also highlighted the existence of a seemingly rogue 
population of apparently older stars, that is, those located well below the main distri- 
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bution (See the HR diagram of low mass stars in Figure 20 left). These have been 
variously interpreted as sources that are coeval with the others but affected by circum- 
stellar material that renders them visible only in scattered light, or sources that truly 
are as old as they appear and therefore offer evidence for large age spreads in star 
forming regions. Although the typical age of ONC members is widely agreed to be 
~l-2 Myr (modulo the systematic effects caused by distance or by adoption of vari- 
ous among plausible sets of theoretical pre-main sequence evolutionary tracks), there is 
still considerable debate regarding the meaning of the apparent luminosity spread. The 
problems illuminated in detail by e.g. McNamai'a, (1976a ) remain. On the one hand, 
there is significant evidence that the error budget for individual stellar luminosities is 
underestimated due to photometric variability, difficulties with extinction corrections, 
unaccounted for multiplicity, etc.; these effects are in addition to known systematic 
problems with the effective temperatures. On the other hand, there are also known pro- 
cesses such as inefficient convection, large amplitude magnetic fields, accretion of new 
material, and perhaps rotation, that pertain to pre-main sequence evolution. Together 
these phenomena suggest caution in any literal interpretation of apparent luminosity 
spreads in HR diagrams as actual age spreads. 

There is, however, evidence from Palla et al. ( 2005 [ 20071) regarding spreads in 
the lithium depletion of young low mass stars in the ONC that seemingly supports the 
range in ages inferred from the luminosity spread in the HR diagram. Ages of 10-30 
Myr are derived for a small fraction of the (notably a moderately to heavily veiled) 
sample. Further, |Jeffries (2007b I used a combination of rotation periods and vsini 
values to infer a distribution of R sin i values and hence a statistical distribution of 
stellar radii R ranging over a factor of 2-3, which they also argue is inconsistent with 
an age dispersion less than 0.3 to 0.5 dex. Clearly more work is needed in order to 
resolve the debate regarding the duration of star formation episodes relative to cluster 
crossing times. 



Stellar Masses and the Initial Mass Function The ONC was an important test case 
regarding ideas of "bi-modal" star formation in which high-mass and low-mass stars 
were suggested as incapable of forming in the same place and at the same time. Empir- 
ically, processes similar to those used to assess stellar ages from HR diagrams were also 



used to assess stellar masses. Evolutionary tracks such as those by Iben & Talbot ( 1966 1 
and Vandenberg & Bell, (1985) showed early on that such a bimodal scenario was not 
applicable to the ONC. Yet it was not until the 1990's that quantitative measurements 
of the initial mass function were claimed for the ONC. 

The first efforts involved use of the surrogate luminosity function, and stellar 
models translated to the empirical plane with some strong assumptions regarding age 



spreads and multiplicity. |Zinnecker & McCaughrean] ( |1991[ see also jZinnecker et al. 
( [T993| )) were the first to apply these techniques to newly available 2 iim survey data 
of the ONC. These efforts were followed by others including [McCaughrean & Stauf- 



fer] ( [T9941 ), |Ali & Depoy| ( [T995l l and |Muench et al.| ( |2000l|2002p that established the 
existence of stars as low in mass as the hydrogen burning hmit (0.073 M0), as well 
as a substantial number of likely brown dwarfs (notably before the confirmation of 
G1229b as the first bona fide brown dwarf!). [Lucas et"aL] ( [2005[ ) have extended these 
arguments to the deuterium burning threshold (0.005 M© or 5 M^jupiter)- Prosser et al. 
([T994) carried out similar analysis using V and / data from HST As an intermediate 
step between one-dimensional luminosity functions and two-dimensional HR diagram 



methods, |Hillenbrand & Carpenter (2000j) applied a statistical technique to H and K- 
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Figure 18. Comparison of IMFs derived for the central 5' surrounding the Trapez- 
ium OB stars in the Orion Nebula. Figure reproduced from Muench et al. (2002i 



band color-magnitude diagrams to determine stellar mass probability functions which 
could be summed to form an initial mass function. 

Hillenbrand (1997) produced the first "forward modeling" method to derive an 



initial mass function in the ONC, making use of the full HR diagram (Figure 17 1 es- 
tablished from optical photometry and spectroscopy, and the evolutionary tracks of 
D'Antona & Mazzitem|([T994]) and |Swenson et al.| ( fl^94l l. [Luhman et al] ( |2000l ) and 
Slesnick et al.| (2004 ) used infrared photometry and spectroscopy to push the investiga- 
tion to lower masses, across the brown dwarf limit. 

All of the above studies are consistent with a mass function that rises in a Salpeter 
like fashion from the highest masses to the sub-solar regime, then begins to flatten 
around 0.5-0.6 Mq with a peak around 0.2-0.3 Mq and then turns over into the brown 
dwarf regime. The exact details depend on the methods, the samples, and the models 



though, as demonstrated by Muench et al. (2002 see Figure 18 1, there is remarkable 
agreement considering possible sources of variance. 



Binarity of T Tauri Stars in the Nebula The Orion Nebula Cluster offers an excel- 
lent opportunity to study the frequency of binarity among young low-mass stars in a 
clustered environment and to determine the shape of the binary separation distribution 
function at an early age. The first major survey for subarcsecond visual PMS binaries 



Pflamm-Altenburg & Kroupa 
than 5Mq. 



( 2006 I argue, however that the ONC suffers a deficit in stars more massive 
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in the ONC was carried out by Prosser et al. ( |1994 1 using the Hubble Space Telescope. 
On their high-resolution optical images they found 35 new visual binaries with sep- 
arations up to 1 arcsec. Padgett et aL] ( [1997| l and Petr et al. ( il998 ) found another 7 
and 4 binaries, respectively. The former study used archived T^-band HST images to 
study an area R > 35" from the Trapezium while the latter used speckle holography to 
reconstruct ground-based near-infra red HK images wit hin the central R < 30". Addi- 
tional visual binaries were found by Simon et al. ( 1999 1 {R <~' 1.5') and Kohler et al. 
(|2006 ) (15 > i? > 5'), both based on adaptive optics K-band observations. The results 
of |Bouy et aLl ( |2008] l illustrate how continued improvements to AO instrumentation will 
add new binary systems to the Orion census. Finally, the COUP X-ray survey identified 
several binaries in the ONC ("Getman et al. 2005a). 

Recently, Reipurth et al., {2001} carried out a major survey of the ONC using 
Ha images obtained with the HST ACIS. In a region extending from 1 to ~ 20' 
(~ 0.1 to 2 pc) surrounding 6^ Ori C they found 78 multiple systems, of which 55 
are new discoveries, with separations less than 1.5". Because of the high stellar den- 
sity in the ONC, it can be statistically determined that 9 of these must be line-of-sight 
associations (see also the analysis of |Simon ( |1997j ) and Bate et aL| ( 1998 1 regarding 
the impact of high stellar density on true binary fraction). When correcting for this, 
a binary fraction of 8.8%ibl.l% is found in the limited range from 67.5 to 675 AU. 



In the same range, the field binary fraction is a factor 1.5 times higher (Duquennoy & 



Mayor|[T991| ), and the binary fraction in loose associations is a factor 2.2 higher than 
in the ONC ( [Reipurth & Zinnecker|1993| . The separation distribution function shows 
unusual structure, with a steep decrease in the number of binaries for separations larger 
than 0.5". Moreover, the ratio of wide to close binaries across the ONC shows a ma- 
jor depression towards the central region, indicating that wide binaries are destroyed as 



they pass through the central potential well, as theoretically expected (e.g. [Kroupa et al. 
[T9^ . 



3.5. Sub-stellar Objects in the Nebula 



Optical and near-IR imaging surveys (Sect. |3.1.| ) have revealed at least 100 sources 
within the central 0.5 pc of the Trapezium that, according to pre-main sequence theory, 
have luminosities at an age of 1 Myr that are consistent with their being substellar 
objects. While similarly deep images of the outer nebula exist, photometric censuses 
have not yet been published. Purely photometric selection techniques are, however, less 
precise at securely identifying brown dwarfs, since their luminosity is a function of age. 
For example, Slesnick et al. ( |2004[ ) reported a number of warmer sub-luminous objects 
that masquerade as lower luminosity brown dwarfs (see Figure 20 right). Whether 
these sources are older or their fluxes contaminated perhaps by scattered light is not yet 
clear. Thus, the identity of an individual young stellai" object (YSO) as substellar rather 
than pre-stellar is best established using measures of that source's effective surface 
temperature and/or surface gravity. The precise boundary at which a YSO can never 
attain sustained nuclear burning comes from theoretical evolutionary models, with a 
working consensus (e.g., Luhman et al.|1998 1 that includes young (<'^ 5 Myr) objects 
cooler than Tgjj <~ 3500 K or a spectral type later than "M6." According to this 
criterion, there are ~ 35 such objects in the most recent version of the Hillenbrand 
(|1997|) wide field catalog. 
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Figure 19. Surface gravity of young Orion brown dwarfs. Left: sharpening of the 
H{1.6 ixm) band spectra compared to field dwarfs with high surface gravities and 
theoretical models. Reproduced from Lucas et al. (2001 1; Right: a comparison of 
the optical Na 1 8 183/8 195 A index for Orion sources, young stars in other regions 
and field dwarfs or giants. Reproduced from Riddick et al. ( 2007[l . 



Surface Gravity Especially over the last five years, the utility of spectroscopic surface 
gravity diagnostics (e.g. Nal, Cal, KI, TiO, VO, CO, steam) has become apparent, es- 
pecially towards the lowest masses, and can in principle be used to distinguish bona fide 
cluster members from both faint red foreground dwarfs and reddened background gi- 
ants. Examples of such surface gravity sensitive observations of Orion Nebula sources 
are given in Figure 19 including Nal (8183/8195A) and steam in the near-IR H band. 
Relevant studies in the optical in clude [Riddick et"aL (2007 see also Figure [T9]). Near- 
see also Figure 19 1, which 



IR studies include Luhman et al. ( 2000 1, Lucas et al, 



(2001 



were the first observations to reveal the sharp "triangle" shape of the H band spectra 
for low surface gravity late M stars, [Slesnick et all] ( |2004[ ), [Slesnick et al.| ( |2005] ), and 
Lucas et aL](|2006). However, there are not yet well established methods at all spectral 



types and the published data on these Orion sources reflect a lack of uniform spectral 
typing. Furthermore, the intermediate-gravity nature of young pre-main sequence stars 
is not readily apparent in all cases. Specifically, many stars with strong evidence from 
both emission lines and infrared excesses for membership are not distinguished as such 
based on surface gravity alone. Establishment of cluster membership may need to rely 
on kinematic association in addition to surface gravity indicators. 



Masses By number, brown dwarfs constitute between 20 and 30% of the total ONC 
membership ( [Muench et al.||2002t [Slesni ck et al.||2004[ [Levine et aL][2006| ), although 
the precise value depends in part upon the assumption of a luminosity and effective 
temperature for a star at the hydrogen burning limit (Levine et al. 20061. By mass 
this is negligible (<'^ 1%) compared to the overall conversion of gas to stars in this 
cloud. The reported substellar IMF alpha (per unit log mass) is ~ 0.6 (Hillenbrand 
& Carpenter 2000; Muench et al.||20()2 i, however, this function is poorly described 
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Figure 20. Left: Central Orion Nebula M dwarf Hertzprung-Russell diagram; 



Right: resulting spectroscopically derived IMF. Reproduced from Slesnick et al. 
( 12^41 [2005] l . 



by a single power-law as seen in Figures 18 and 20 (right). Additionally, the near-IR 
luminosity function of the central ONC displays structure at the faint end that is greater 
than expected purely from field star contamination but is not consistent with a simple 
declining power-law IMF given current PMS theory ( [Muench et al.|20d2 ). 



Interesting Very Low Mass Sources The object 2MASS J05352 184-0546085 was 
found to be a unique brown dwarf-brown dwarf double line eclipsing binary system 
by Stassun et al. (2006a). The combination of photometric and spectroscopic monitor- 
ing allowed individual masses (~ 0.04 M0) and radii (~ 0.6 Rq) to be determined. 
Interestingly, the higher mass component has a cooler surface temperature than the 
lower mass component (Stassun et al. '2007a', '2008*); more recent evidence of strong 
surface activity on the primary (| Reiners et al.,2007j ) indicates this temperature reversal 
could be explained by the inhibition of convection by strong magnetic fields (D'Antona 
et al.||2000 ). The sustained efforts of deep monitoring projects in Orion are beginning 
to uncover additional low mass double line eclipsing systems, which can provide im- 
proved constraints on very low mass pre-main sequence evolutionary theory. These 
include JW 380 ( jlrwin et al7||2007| M ~ 0.2 Mq) and Parenago 1802 ( |Cargile et al 
2008| M ~ 0.4 Mq). 



3.6. Circumstellar Matter 

With a rich stellar population spanning the mass spectrum all the way from 0-type 
massive stars to late-M-type brown dwarfs, the ONC is an obvious arena for investiga- 
tions of circumstellar disks. Known for decades as a rich collection of emission line 
objects (e.g., Haro|1953 Parsamian & Chavira|1982 1, the ONC was also an early target 
of infrared investigations. In particular, the single channel photometer measurements 
of IRydgren & Vrba| ([T9841), [Breger efal] ([T98T]), |Smith| (|T976ll, |McNamara| ( |1976a| ). 



PenstonI p973] ); [Penston et aL| ( [T975| ), |Ney et al.| ( [1973) , and |Lee] ( [1958) established 
that many ONC stars had excess emission at near- to mid-infrared wavelengths. The 
frequency of circumstellar disks around brown dwarfs in the Trapezium appear to be a 
smooth continuation of this property for higher mass objects ( [Hillenbrand et al.|1998 
Lada et al.|2000t [Muench et al.|2001t|Lada et al.|2004) . 
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The advent of infrared arrays led to the studies of McCaughrean & Stauffer|(|1994 ); 
McCaughrean & O'Delli (fT996b, |Ali & Depoy' fl995), 'Jones et ai; (1994), 'Hillenj- 
brand et al. (1998), Simon eraL|([1999[),|Hillenbrand & Carpenter (2000) , Lada et al. 
( |2000[ 12004) . IMu'ench et al.| ( |2001| ), [Lucas & Roche] ( |2000| ), which provided a cen- 
sus of JHKL excesses, an assessment of the prevalence of circumstellar disks with 
stellar mass, and correlations with stellar environment within the ONC. Most recently, 
Spitzer data has improved the detail ed knowledge of thes e dis ks at mid-infrared (3.5 — 
24 /im) wavelengths with papers by |Rebull et al.|(|200 6') and' Cieza & B aliber'('2007[) 
adding to ground-based mid-IR work, e.g., Stassun et al. (2001 1, in investigating disk 
versus stellar rotation paradigms. Additional mid-infrared work focussing on proplyds 
is that of Smith et al. (200 5 [). Towards longer wavelengths, millimeter investigations 
including those of ^Eisner & CarpenTer| ( |2006| ), [Eisner et al.| ( |2008] ) and [Williams et al. 



( 2005 1 have together established the prevalence of minimum-mass solar nebula disks 



in the ONC, and the existence of several much more massive systems. 



Disk Accretion Studies of the accretion signatures typically associated with such 
disks are difficult in the ONC due to the bright continuum and emission-line backdrop 
of the Nebula. However, studies of the Call "infrared" triplet ( Samuel[[l993 ; Hillen 
brand et al. 1998), which is the only strong optical line found in accretion disks but 
not in HII regions, avoids issues with strong and variable background in more typically 



observed lines such as Ha. Using high dispersion spectroscopy, Sicilia-Aguilar et al. 



(2005 1 and Filresz et al. (2008 1 were able to separate nebular and stellar emission at Ha 



and, focussing on the line wings, claim to find many sources (including those lacking 
thermal excess) that indicate accretion rather than winds. Finally, studies of the accre- 



tion continuum at UV wavelengths have included those of Robberto et al. ( 2004 1 in the 
inner ONC using HST and Reb uU et al.| ( |2000l ) in the ONC "flanking fields" from the 
ground. In summary, 80 — 100% of the stars in the inner ONC are suspected accretors. 
A median accretion rate of 10^^ M^/yr has been derived, with some evidence for a 
dependence of accretion on stellar mass. Deep spectroscopic studies of such accretion 
signatures are lacking for the brown dwarf population. 



X-ray Clues to Disk Evolution The high X-ray intensity and hard spectra found for 
ONC stars indicate that the ionization of disk gases by stellar X-rays dominates ioniza- 
tion by cosmic rays or other sources by a large factor (~ 10^ for 1 M0 stars). COUP 
provides two lines of direct evidence for establishing where the disk irradiation by X- 
rays originates. First, Kastner et al.[ ([2005 [) find that the X-ray absorbing column of 



the COUP stars surrounded by dusty "proplyds" imaged with HST increases with disk 
inclination. The soft X-rays being absorbed must ionize disk gas. Second, the COUP 
results(and Chandra studies of other YSOs; e.g. Imanishi et al.|200 l) include the de- 
tection of the fluorescent emission line from cold iron atoms at 6.4 keV, which is seen 
next to the hot plasma line around 6.7 keV (Figure [21] Tsujimoto et al. 2005[ ). The 
equivalent widths of the 6.4 keV line observed in the COUP sample are compatible 
with the fluorescence originating in a centrally illuminated disk observed face-on, and 
cannot be attributable to the fluorescence by interstellar or circumstellar matter along 
the line of sight. Recent theoretical studies have predicted that the behavior of proto- 
planetary disks and the processes of planet formation will be significantly altered if the 
disks are slightly ionized (see reviews by |Feigel son[[2005[ [Glassgold et al.[[2005D. For 
example, the planets would form in a turbulent and lumpy disk rather than a smooth 
disk, which may prevent Earth-like planets from rapidly migrating through the disk 
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towards the young star. Thus the X-rays from young stars may have important imphca- 
tions for the formation of planets around these stars. 
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Figure 21. The signature of disk irradiation by X-rays from powerful flares of 7 
COUP sources. The presence of a fluorescent emission line from cold iron atoms at 
6.4 keV is seen next to the hot plasma line around 6.7 keV (both marked in the figure 
with arrows). Reproduced from Tsujimoto et al. (2005 1. 



3.7. Kinematics and Cluster Structure 



Proper Motion Studies Proper motion studies in the literature are greatly muddled by 
the quality of and by systematics evident in the observational data. This has yielded 
conflicting results about the relative motions of stars and whether they display system- 



atic expansion or contraction. Authors including, Parenago ( 1954a I, Strand ( 1958 ) and 



Fallon et aL] ( |1977| ), have all produced evidence for the expansion or contraction of 
the cluster, with evidence offered by [Vasilevskis ( 1962 1971 1 and Allen et al. ( 1974 1 
indicate that these claims are due to observational error. Two classic proper motions 
studies are |van Altena et al.| ( [T988| ) and [Jones & Walker] ( [T988p . The [van Altena et al. 
proper motion study used plates from three observatories over a 77 year period. The 
internal error estimates are very small for the high mass members {V < 12.5) and they 
derived a velocity dispersion of 0.7 milliarcsec/yr (1.3 km s~^ at a distance of 400 pc) 
for 48 members. [Jones & Wal ker| (|1988]l used 47 deeper red plates taken over a 20 
year timescale to derive proper motions and membership probabilities for nearly 1000 
stars near the Nebula. Focusing on somewhat lower mass objects, they find a velocity 
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dispersion of 2.5 km s~^, which they point out is similar to the clump to clump veloc- 
ity dispersion of the gas. They also confirmed van Altena et al.'s finding of a smaller 
velocity dispersion for the bright stars. This suggested to Jones & Walker a relation- 
ship between velocity dispersion and mass that could also constrain the dynamical state 
of the cluster. Reviewing these results, Hillenbrand & Hartmann ( 1998^ find that any 
variation in velocity dispersion with mass appears to be too small to be consistent with 
the equipartition of energy occurring during significant dynamical evolution. The long 
baseline and high resolution of data in the HST archive should provide good astrometry 
for updating these proper motion results (e.g., O'Dell et al.^lOOS) . 
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Figure 22. Distribution of (heliocentric) radial velocity as a func tion of declina- 
tion for stars (dots) and gas (^■^CO as the blue background, from Bally et al.|l987| l 
in and around the Orion Nebula. Reproduced from Furesz et al. (2008 ). Each panel 
represents a different range of right ascension; the color and sizes of the dots cor- 
respond to the net precision of their radial velocity measurements (larger or black 
dots being the best quality); the last panel compares the stellar and gas distribution 
integrated across the cloud. 



Radial Velocity Studies The difficulty of acquiring radial velocities for a large number 
of young stars has been recently overcome using new multi-fiber echelle spectrographs. 
While earlier works such as Johnson] ( |1 965 j) or|Walker|(|1983|) w ere able to observe tens 



of stars, the recent surveys published by Sicilia-Aguilar et al. (2005 1 and Fiiresz et al. 



(|2008) include a total of about 1300 stars around the Orion Nebula (and the Northern 
Orion A Molecular Cloud). They find a bulk velocity dispersion of ~ 2.3 — 3 km s~^ 
for all the members. This dispersion should be interpreted carefully considering the 
existence of a large ~ 5 km s"^ gradient across the Nebula, which is most evident in 



the top right panel of Figure 22 Indeed across the Northern Orion A cloud, Fiiresz et 
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al. find that the radial velocities of the stars appear to correlate strongly with the radial 
velocity of the molecular gas cloud. 

Structure The Orion Nebula Cluster is elongated parallel to the Orion A Molecular 
Cloud with a centroid from elliptical star count fitting falling just North of 9^ Ori C 
( [Hillenbrand & Hartmann 1998) . No substructure is evident in the optically revealed 
stars ( Hillenbrand & Hartmann 1998} [Scally & Clarke 2002 ), which has been inter- 
preted to suggest that the cluster may have formed from a fairly large number of small 
clusterings that quickly dispersed. Some evidence for subclustering at the < 0.1 pc 



scale appears when one studies the youngest, most embedded members (Lada et al. 
[20001 [20041 [Grosso et ar[[2005] ). Hillenbrand & Hartmann also found that the cluster 
displays a radial profile that is well fit by a King profile, ha ying a core radius of 0. 1 
0.2 pc and a central stellar density exceeding 10^ stars • pc~^ ( McCaughrean & Stauffer 
[1994; Hillenbrand & Hartmann 1998). It has been argued that this mass distribution is 
primordial as it seems unlikely that the cluster achieved a King like radial profile as the 
result of dynamics given its youth. A similar, primordial explanation is given regarding 
evidence for mass segregation of the highest mass stars to the cluster's core and a skew 
of lower mass stars to the outer parts of the cluster ( |Hillenbrand & Carpenter[[2000 1. 
Investigations into the spatial distributions of stars and brown dwarfs have been made 
( jKumar & Schmejaf2007[ ), but that field of study was very small and perhaps more in- 
teresting results await much wider field studies that cover more than the central 0.5 pc 
and include radial velocity studies of the faintest members. 



4. The Most Massive Stars in the Orion Nebula 

The Orion Nebula is an HII region created (primarily) by the ionizing photons of the 
central O star, 9^ Ori C. There are a total of 20 optically revealed O or B type pri- 
maries, which are listed in Table [4[ this tabulation is a merger of Orion 1 subgroup "d" 
association members from Brown (1996 ) and the Hillenbrand (1997 ) ONC census. 

The collective Bayer designation for the Orion Nebula is 9 Orionis. The nomen- 
clature of members of the inner Orion Nebula break down further into 9^ and 6*^ desig- 
nations, where 9^ Ori is the Bayer designation for the famous "Trapezium" asterism. A 



finder chart for the Trapezium is given in Figure 23 The principal components of the 
Trapezium are further labeled A, B, C, D, E, F; e.g., 6^ Ori C; subcomponents of these 
'A,B,C..." designations are marked with subscripts, e.g., 9^ Ori Ci. Five of these six 
primaries in 9^ Ori are O or B type stars. The 9"^ Ori designation corresponds to a string 
of three OB stars, e.g. the A,B,C components, lying near the bright bar; they are eas- 
ily identified in Fig[7[just southeast of the Trapezium. Additional cross-references and 
observed positions for the unresolved primaries of these systems are given in Table [4[ 



The identifiers from the Brun ( 1935 1 catalog especially useful for interpreting older 
texts. For completeness the star names 41 and 43 Ori are the Flamsteed designations 
for 9^ Ori and 9^ Ori A sources, respectively. Finally, the Trapezium is also frequently 



found listed by its catalog entry in the [Aitken[ ( [1932[ ), ADS 4186. 
4.1. Basic Properties 



Simon-Dfaz et al.|( 2006| l derive MK spectral types for 5 of the OB members in the cen- 



tral nebula and provide stellar properties (R, T, M, L, g, v sin i) and extensive analysis 
of oxygen abundances through comparison to synthetic stellar models and template OB 
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Table 4. Massive (O or B spectral type) stars of the Orion Nebula: Names, Cross- 
References and Positions. Positions of the primary star are given in the mean epoch 
of their observation except where marked ("*")■ Reference catalogs: "2MASS" 



Skrutskie et al.|(l2006l l; "C OUP" [Getman et al.| ( p005cl l; "HIP" (Hipparcos) [Pe^ 
man et al.| ( |1997| ; "MLLA " Muenc h et al.|p002| l; "PPM" |Roeser & Bastian| ( |1988| l; 
"TYC" (Tycho-2) |H0get al., ( ,2000i ) 



Designations Positional Information 



Parenago Common Draper Brun Reference Catalog Equatorial (ICRS) Epoch 

(1954) (HD) (1935) Name # RA DEC (yr) 



1539 






328 


2MASS 05343999-0510070 83.666636 


-5.168633 2000.9 


1605 


V 327 Ori 


36917 


388 


TYC 


4774-809-1 


83.69575750 


-5.57072639 1991.7 


1660 




36939 


442 


TYC 


4744-823-1 


83.73038972 


-5.50614139 1991.7 


1744 




36981 


502 


TYC 


4744-915-1 


83.77582389 


-5.20442611 1991.7 


1772 


LPOri 


36982 


530 


TYC 


4774-849-1 


83.79098528 


-5.46478444 1991.7 


1865 


6*1 Ori A 


37020 


587 


HIP 


26220 


83.81592896 


-5.38731536 1991.25 


1863 


el Ori B 


37021 


595 


COUP 


778 


83.817270 


-5.385220 2003.04 


1891 


6*1 Ori C 


37022 


598 


TYC 


4774-931-1 


83.81860444 


-5.38969611 1991.8 


1889 


6*1 Ori D 


37023 


612 


HIP 


26224 


83.82166016 


-5.38768076 1991.25 


1892 


6*1 Ori F 




603 


MLLA 


388 


83.819625 


-5.390222 2000.3 


1956 


V 1230 Ori 




655 


PPM 


702316 


83.83632500 


-5.36236111 2000(*) 


1993 


Ori A 


37041 


682 


TYC 


4774-933-1 


83.84540917 


-5.41606333 1991.7 


2031 


6*2 Ori B 


37042 


714 


TYC 


4774-934-1 


83.86000444 


-5.41687306 1991.7 


2085 


Ori C 


37062 


760 


PPM 


188231 


83.88098334 


-5.42122222 2000(*) 


2074 


NUOri 


37061 


747 


TYC 


4774-906-1 


83.88068417 


-5.26738778 1991.7 


2284 




37114 


920 


TYC 


4774-867-1 


83.99391639 


-5.37537417 1991.7 


2271 




37115 


907 


TYC 


4778-1369-1 


83.97533639 


-5.62839861 1991.6 


2366 




37150 


980 


TYC 


4778-1378-1 


84.06261139 


-5.64792000 1991.6 


2387 




37174 


992 


TYC 


4774-855-1 


84.11326917 


-5.40870417 1991.7 


2425 






1018 


TYC 


4774-873-1 


84.15876833 


-5.47638250 1991.7 



Stars. An additional source of basic stellar properties can come from the multi-spectral 
flux ratios of the various binary components. For example, Kraus et al. ( 2007 1 used flux 
ratios from V to K bands for the 9^ Ori Ci^2 binary to constrain the components to have 
spectral types of ~05.5 and 09.5. Papers surveying rotational and radial velocities (to 
the extent permitted by the high multiplicity) of these objects include |Abt et al.| ( |l970| l, 
Abtet all ([1991) and |Wolff et aT1 ( |200gi ). 



X-rays from the OB stars were thought to be generated in a myriad of tiny shocks 
in their powerful winds. But this model predicts that the X-ray spectrum will be soft and 
the emission will be constant. The Chandra Orion Ultradeep Project (COUP) shows 
that the spectra of the Trapezium OB stars often have a hard component and can exhibit 
rapid flares ( [Stelzer eTa l. 2005 ). Of the 10 unobscured COUP sources earlier than 
B4, only three (9^ Ori D, NU Ori, and possibly 9^ Ori B) show the expected signature 
of many small wind shocks, while most show flares and/or hard spectral components. 
This suggests that the winds of these OB stars are, at least in part, trapped by magnetic 
fields, resulting in large scale shocks and production of hard X-ray emission ( jBabel & 



Montmerle 19971. Figure 24 compares X-ray emission from the steady-wind source 
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Figure 23. Finder chart for the Trapezium and the inner Orion Nebula. North is 
up and East is to the left. 



9^ Ori D with that of the magnetically active 6^ Ori A. COUP also confirms that the 
X-ray emission from intermediate-mass stars with spectral types B5-A9 is attributable 
to lower mass companions rather than the intermediate-mass star itself. Additional X- 



ray studies of the OB stars in the Orion Nebula include Caillault et al. 



etaL]([T995l ) IStahl et al.|dT996l),|Schulz et al.| ( [20001 ), ISchulz et al.|^200T] l, |Schulz et al. 
(|2003aD andlSchulz et al. (|2003b|l. 



(19941, Geier 



4.2. Kinematics 

Compact stellar groups like the Trapezium are expected to be intrinsically unstable 
( Ambartsumia n|1935 1. Therefore, the kinematics of such a stellar grouping could dis- 
play systematic expansion or contraction, display evidence for disintegrating multiple 
systems or provide an origin for runaway stars ( Poveda et al.|1967 ). Regarding signa- 
tures of disintegrating systems, [van Altena et al. '{ 1988 1 reported that the two O stars in 
the central Orion Nebula, 9^ Ori C & Ori A, were observed to have proper motions 
relative to other bright {V < 12) stars that were large enough to carry them out of the 
center of the Nebula in less than 1 Myr. However, the relative proper motions for the 
Trapezium stars derived by ,Allen et al., {1974) and updated in Allen et al., ( ,2004) are 
in general very small, which calls into question whether 9^ Ori C is actually be ejected 
relative to the other Trapezium OB stars. 

The dynamics of two more systems related to the Orion Nebula have garnered 
much more detailed attention. The first is the existence of two high velocity OB stars 



whose origins appear to coincide near the Orion Nebula. Blaauw & Morgan ( 1953 ) de- 
scribed the high velocity and apparent space motion of AE Aurigae away from Orion, 
while Blaa uw & Morgan] ( | 1954 ) made this interesting connection for the star fi Colum- 
bae as well. While each is over 25° from the Nebula, they have proper motion vectors 
corresponding to space motions > 100 km s^^ and Blaauw & Morgan (1954) esti- 
mated that they both originated in a dissolution event occurring near to the current 
Nebula about 2.7 Myr ago. Subsequent analysis included |Blaauw| ( [ 1 96 1 1 ), while Gies 
& Bolton[ ( |1986| ) proposed adding the O star i Ori, which is just South of the Nebula, 
into the dissolution. [Hoogerwerf et al. ( 2000[ ) and Hoogerwerf et al.| (|2001|) performed 
the integration of the orbits including improved Hipparcos data, radial velocities and a 
Galactic potential field; their results further constrained the dissolution event to have 
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0M517.J-0S2J16 




lU 200 



Figure 24. Top: lightcurve and spectrum for 9^ Ori D with the soft constant emis- 
sion, as a signature of many small wind shocks. Bottom: in contrast to 9^ Ori D the 
lightcurve and spectrum of 9^ Ori A show hard flaring emission possibly due to the 
confinement of the wind by a strong stellar magnetic field. In the lightcurves, the 
black line indicates the total (0.5 — 8.0) keV energy band, red line - soft (0.5 — 2.0) 
keV band, and blue - hard (2.0 — 8.0) ke V band. Chandra X -ray spectra and 
lightcurves from the COUP. Reproduced from Stelzer et al. (2005[l. 



occurred approximately 2.5 Myr ago and could have been located nearer the center of 
the cluster even than the location of t Ori today. Numerical modeling of this event has 
been performed most recently by Gualandris et al.| (2004l. 

The second apparently disintegrating high mass system concerns the relationship 
between the B-N object and other high mass protostars or O stars in the central cluster. 
The relatively high proper motion of the B-N object relative to the cluster was derived 
first by'Plambeck et al. (1995); monitoring of the proper motion of this object using 
published ( ^Rodriguez et al.,2005, Gomez et al.,2005 ) and unpublished data have led to 
two hypotheses on the origin of its motion. Figure 25 compares proper motion diagrams 
tracing the B-N object back to an origin in two different formation scenarios. (Tan| 
( |2004) argued this is a recoil motion from the B-N object being ejected from 9^ Ori C 
~ 4000 years ago, while Gomez et al. (2005 1 argue for a dissolution of the B-N object 
from the protostellar sources "i" and "n" a few hundred years ago. Either scenario 
claims to be consistent with the production of the explosive outflow in the K-L region 
500-1000 years ago. (see O'Dell et al., Part II). While this outflow is a consequence 
of the dissolution of the B-N, "i" and "n" system in the Gomez et al. scenario, the 
hypothesis of Tan is that the B-N triggered the outflow by passing extremely close to 
source "i." 

Recent high (< 0.1") optical/near-IR observations have allowed the detailed or- 
bits of the many multiple high mass systems (see next subsection) to be constrained. 
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Figure 25. Proper motion histories for observations relating to the origin of the 
B-N object as a runaway B star. Left: t he d i sintegr ation of the B-N object and 
6^ On C 4,000 years ago as described by 'Tan (2004i; Right: the disintegration of 
the protostellar triple B-N, source "i" and source "n" 500 years ago as described by 
(Gomez et aL|(201)5] ). 



These include the recent works by'Schertl et al. (2003 1, Close et al. ( |2003| l, Kraus et al. 
( [2007^ , and Patience et al., ( ,2008j . There is, however, an almost complete lack of such 
high resolution imaging or monitoring of high mass objects outside of the immediate 
Trapezium core; such observations could be very valuable for better understanding the 
kinematic history of the high mass stars, their interactions, and the dissolution of higher 
order systems. 



4.3. Multiplicity 

The massive members of the ONC have been the target of several studies concerning 
the multiplicity of these stars. Light curves revealed a number of eclipsing systems 
(e.g. [Parenago & Kukarkin|[T947| [Hall & Garri"sonl[T969l |Lohsen|[T975| |Wolf||1994|l, 
which have been further analyzed to identify additional companions (e.g Vasileiskii & 
Vitrichenko|200q 



Searches for spectroscopic companions were presented by,Abt et al. 



( 1991 1 and Morrell & Levato ( 1991 1, who found spectroscopic binary frequencies of 
20% - 30%. Studies that detected visual companions to massive ONC members include 
Padgett et al.] ( fT997l ), |Petr et al.| ( fT998] ), and |Simon et"aL] ( fT999] ). [Preibisch et aLl ( [T999l l 
performed a systematic survey for multiple systems among 13 bright ONC members 
of spectral type O or B with the technique of near-infrared bispectrum speckle inter- 
ferometry, and complemented the results with information about known spectroscopic 
companions. In the speckle images, which have a resolution of 0.075", eight visual 
companions were found (see Figure [26] ). Stellar masses of the companions were es- 
timated from the observed near-infrared flux ratios. The properties of these multiple 
systems are summarized in Table |5] 

The results of this multiplicity survey allowed general conclusions about the mul- 
tiplicity of the massive stars in the ONC. The mean number of known companions in 



45 



M jUiplic:i1y in tti? Orion Trapezium 

W^lgiiLeCdL \a99: 2O07) 




Figure 26. The multiplicity of the stars in the Trapezium systems as revealed by 
infrared interferometry. 



the sample of the 13 observed OB stars is 15/13 = 1.2 per primary star. Since it is 
highly likely that there are even more, still undetected companions, this is clearly a 
strict lower limit for the true multiplicity. The multiplicity of these stars is thus con- 
siderably (at least about 2 times) higher than that among low-mass primary stars in the 
general field as well as in the ONC. Multiplicity seems to be a function of spectral 
type: based on the literature summarized above stars of spectral type B3 or later have, 
on average, 0.6 known companions per primary, while the stars of spectral type O to 
B2 have at least 1.5 companions per primary. This comparison clearly suggests a mass- 
dependence of the multiplicity, which is in qualitative agreement with the generally 
observed trend of increasing multiplicity with increasing stellar mass seen among low- 
and intermediate-mass stars. While most low-mass stars may be single (|Lada'2006), the 
data suggest that essentially all high-mass stars are members of higher-order multiple 
systems (triples, quadruples). 

These results point towards a fundamental difference between high- and low-mass 
stars and support the assumption that the formation of massive stars is not simply a 
scaled-up version of low-mass star formation Zinnecker & Yorke (2007). Theoretical 
scenarios for the formation of these massive stars, including stellar interactions, the cap- 



ture of companions and perhaps even stellar mergers scenarios are discussed in Bonnell 



etaL](|2004t;|Bo nnell & Bate|([2005]) [Moeckel & Bally| ( |2006l |2007a|b| ), [Davies et al. 
p006), and|Bally & Zinnecker p005]). 
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Table 5. Massive stars of the Orion Nebula: multiplicity. Adapted from the sur- 
vey of Pre ibisch et al.| ( |T999| l. The Table lists the name, spectral type, and mass of 
the primary component, and then the nature (visual or spectroscopic companion), 
separation, and mass ratio of the known companions. If there is a constraint on the 
existence of a binary companion, then that limit and its reference are given. Refer- 
ences (Column Ref) correspond by number to the following s ources: 1: ,PreibiscTi| 
"etaL]([T999|; 2 : Weigeh eTaL](fT999]l; 3: iPetr et al.|(|T998); 4: |Bossi et al.|(|1989|l; 



5: ISimon et aT](11999j; 6:lAbt et al.|(|199l| ); 7 
P007»; 9|Vitrichenko & Plachiiida|(200Tf; 10: 



ILevato & Abtl(|1976| ); 8: |Kraus et al.| 
al.|(|2006D 



Kohler et ; 



Parenago 


Other 


Mp 


Mt 
r A/f 1 


Companion 


P 

r A T n 


M/Mp 


Ref 


1605 


V 327 Ori 


3.5 


7 


-2 (spec) 




~ 0.9-1.0 


7 


1660 










< 60 


— 


10 


1744 










< 60 




10 


1772 


LP Ori 








< 100 




1 


1865 


Ori A 


14 


21 


-2 (vis) 
-3 (spec) 


100 
1 


~ 0.25 
0.19-0.33 


2,3 
4,9 


1863 


Ori B 


7 


12 


-2 (vis) 
-3 (vis) 
-4 (vis) 
-5 (spec) 


430 
460 
260 
0.13 


~ 0.22(<0.71) 
~ 0.10(<0.50) 
~ 0.03(<0.29) 
~ 0.39 


1,2 
1,2 
1,5 
6 


1891 


Ori C 


34 


50 


-2 (vis) 


18 


~ 0.45 


2,8 


1889 


Ori D 


16 


> 16 


-2? (vis) 


8 




8 


1993 


Ori A 


25 


40 


-2 (vis) 
-3 (spec) 


173 
0.47 


~ 0.28(<0.32) 
~ 0.35 


1 
6 


2031 


e"^ OriB 








< 100 




1 


2074 


NUOri 


14 


17 


-2 (vis) 
-3 (spec) 


214 
0.35 


~ 0.07(<0.28) 
~ 0.2 


1 
6 


2284 










< 60 




10 


2271 


HD 37115 


5 


7 


-2 (vis) 


400 


~ 0.29(<0.96) 


1 


2366 










< 100 




1 


2425 




4 


5 


-2 (vis) 


388 


~ 0.04(<0.35) 


1 



In the near future, long-baseline interferometric surveys of the ONC OB stars will 
allow resolution of companions as close as 0.001" (0.45 AU). This will largely fill the 
still existing gap between currently detectable visual and spectroscopic companions and 
provide much more complete information about the multiplicity of these stars. 



4.4. Notes on Individual Trapezium (O^Ori) Stars 

Ori A; Parenago 1865 (V1016 Ori; HR 1893; HD 37020) Th e bright {V = 6.7) 
westernmost member of the Trapezium, of spectral type B0.5V ( Simon-Diaz et al^ 
2006| ), was discovered as recently as 1974 to be an eclipsing binary (Lohsen 1975| l 
with a period of 65.43 days and deep primary minima (AV ~1); no secondary min- 



ima have been seen with certainty. Early attempts to model the system include Bossi 
et al.| ( |1989| ), |Vitrichenko et al.| ( [1^981 ) and |Vitrichenko| ( |1999| ). Optical and ultraviolet 
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spectra are analyzed in 



Vitrichenko & Klochkova 



richenko & Plachinda (2001 1. Stickland & Lloyd ( |1999[ ) review and analyze existing 



(200U), '^^chenko| ( |200T] | and[W 



data up to that time, deriving masses of 12 and 3 Mq for the eclipsing pair. ,BondaF| 
et al. (20001 analyzed new light curves, and Vitrichenko & Plachinda] ( [2001 1 suggested 
that the components have masses of 21 Mq and 3.9 Mq. A primary mass of 14 Mq 
and a mass ratio q 0.2 seems a reasonable summary of rec ent finding s . A th ird star 



Petr et al. 



(19981; in the 



(A2) at 0.2" separation has been found next to 9^ Ori Ai 3 by 
near-IR its flux ratio to the eclipsing primary suggests a reddened A or F type pre-main 
sequence star. Radio continuum observations of V 1016 Ori at 2 and 6 cm (|Church 



well et al.|[T987| [Felh et al.|[T989l [T99T| |1993b| ) reveal a variable non-thermal source 



now confidently associated with this visual secondary ( Petr et aL]|1998 ; [Petr-Gotzens 



[& Massi„2007 ); the system is unresolved in Chandra data. Various explanations for 
this radio emission have been proposed, including wind-wind collisions. Orbital mo- 
tion between A2 and the eclipsing primary Ai^3 has b een observed through multi-epoch 
observations over a 7 year period ( [Close et al [2003 ; Schertl et al.|2003 



Ori B; Parenago 1863 (BM Ori; HR 1894; HD 37021) This source consists of at 
least 5 stars but possibly as many as 7. The long-known visual companion, at a separa- 
tion of about 1" from the primary Ori Bi, is clearly resolved into a close (0.117") bi- 
nary (e'lQriBas). At least one member of the B2,3 system is a proplyd, source identifier 
160-307 (note that lO'Dell & Wong| ( [T996 ) gave'this source the identifier 161- 307 and 



considered it stellar); both appear extended in the bispectrum band images of Schertl 
etal.| ( |2003| l; one of them drives a microjet (HH 508) (Ball y et al.|2000|); in the thermal 



infrared, the 6^ Ori B2,3 system dominates the Bi primary (Lada et al. 



2004 



Smith et al. 



2005 1. This system is also the dominant X-ray source and exhibits all the features of 



typical solar mass T Tauri stars including multiple flares ( Stelzer et al.|20()5 1. Another 



faint v isual companion, 6^ Ori B4, is located 0.578" from the primary star (Simon et al. 



1999 1. The primary, 6*^ Ori Bi, is itself an eclipsing binary with a period of 6.46 days; 



early observations are discussed by |Hartwig| ( |1921a|b| l, [Parenago & Kukarkm ( |1947| l 



and Schneller ( 1948). A near-IR light curve is shown in Figure 27 (Sect. 5.1.1. This 



system is also a double-lined spectroscopic binary, so masses and radii of the two com- 



ponents can be derived. The primary is a B3 star. Hall & Garrison ( 1969 1 obtained a 
detailed light curve, demonstrating that the primary eclipse is total with a duration of 9 
hours, and finding a shallow second minimum. Spectra obtained by Popper & Plavec] 
( 1976 1 revealed that the secondary component is a late A-type star. A primary mass of 
~ 7Mq and a mass ratio g ~ 0.4 is adopted for this system. New light curves obtained 



by |Wolf| ( [T994^ were analyzed by |Vasileiskii & Vitrichenko] ( |2000[ ), who suggested that 
there is a third (B type) star in the primary 6^ Ori Bi^s system that is not involved in 
the eclipse, but affects the spectra observed. Vitrichenko et al. ( 2006] ! found a fourth. 



late-type component in the primary system based on a radial velocity anomaly (Popper 



& Plavec|1976[ 


Vitrichenko & Klochkova|2004 


Malov 


(2006 


1 detected high-velocity outflowing 



Ultraviolet spectra by Vitrichenko & 



provide the best-resolution observations with a wide field astrometric reference frame. 

9^ Ori C; Parenago 1891 (HR 1895; HD 37022) The most massive star in the ONC, 
is an extremely close visua l binary system with an initial disc overy o f the companion 



at a separation of 0.033" by ]WeigeU et al.](]1999|. ]Schertl et al.](]2003]) subsequentl y re- 
ported the detection of orbital motion. Kraus et al. '{ 2007 1 and Patience et al. ( ]2008j ) re- 
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port recent multi-epoch observations with visual and near-infrared bispectrum speckle 
and near-infrared long-baseline interferometry. The current data trace the orbital mo- 
tion of t he companion over a more than 10-year period and cover a significant part of 
its orbit. Kraus et al. (20071 derived a highly eccentric (e ~ 0.91) orbit with a period 
of 10.9 yrs and a total mass of ~ 50AIq with the primary being a 34 Mq 05.5 star 
and the companion a 15.5 Mq 09.5 star. The addition of 6 epochs over a baseline of a 
year by [Patience et al. (2008 1 suggest a longer period 26yr), a less eccentric orbit 
(e ~ 0.16) and a somewhat lower total mass (~ 4OM0). 

Another source projected very close to 9^ Ori C is a mid-infrared source found 
~ 2" West of the O star. A ring-like structure around this object is evident in thermal 
and mid-IR images ( Lada et al.|2004 Smith et al. 2005 1 and it is a VLA radio source. 
Source identifiers for this object includ e VLA- 16 (,Felli et al.||1993b[ ), SC3 ( [Hayward 
et al.||1994|l, the L band source #268 (ILada et al. 1200411, and the mid-IR identifiers 



MAX 106 ( |Robberto et aH2 005 ) and 163-323 ( [Smith et al. '2005 ). The spectral energy 
distribution of VLA-16 was modeled as an irradiated proplyd by Robberto et al. ( [2002 j ), 
though it is not entirely clear where the source lies along the line of sight. 

Additional radial velocity components that do not turn out to correspond to this 
binary pair were reported by Vitrichenko] ( 2002 1. Understanding such radial velocity 
variations in 6^ Ori C is complicated by the source's strong magnetic activity, its stellar 
winds and a 15 day periodicity apparent in radio and radial velocity monitoring. These 



winds and strong magnetic field have been investigated b y[Walbom & Nichols ( 1994 



( [Simon-Dfaz et al.|2006| l, [Wade et al.[ ( [2006[ ), and most recently [Stahl et al.[ ( [2008[ ). 



Stahl et al.[ ([T996]l, [Ponati et al.[ ([2002|), [Gagne et al.[ ( [2005[ ) , ([Stelzer eTal 



2005 



9^ Ori D; Parenago 1889 (HR 1896; HD 37023) The second most massive primary 
in the Trapezium has a B0.5V type and a mass of 16 Mq (Simon-Diaz et al. [2006| . 



The source display s a constant soft X-ray spectrum (lacking flares or hybrid spectra) 
([Stelzer et al."2005), suggesting it is neither magnetic nor contains a low mass [q < 0.2) 
companion. While Kraus et al. (2007) found evidence in the UV plane for a companion 
at 8 AU separation, having 10% of the luminosity of the primary, they also point out this 
could be the result of a disk's inner edge. 9^ Ori D is surrounded on 2 sides by the Ney 



& Allen nebula ( [Ney & Allen,, 1969^ |McCaughrean & Gezari„1990) , a 10-20 micron 
diffuse structure whose origin is unclear. Robberto et al. ( 2005 1 present evidence that 
this dust arc cannot originate simply from the ionizing photons of 9^ Ori C but requires 
a wind from 9^ Ori D and a substantial reservoir of dust, which they infer to arise in 
a disk. On the other hand. Smith et al. ( 2005 1 propose that the dust shell is the result 
of the winds of 9^ Ori D interacting with the PDR as it plows down into the remnant 



molecular cloud ( |O'Dell|200l] ). 



9^ Ori E; Parenago 1864 The fifth brightest member of the Trapezium is 9^ Ori E, a 
very strong X-ray source, the second-strongest in the central cluster after 9^ Ori C (Ku 
et al.|1982 i. FelU et al. ( 1993a|b 1 found the star to be a bright and variable non-thermal 
radio co ntinuum source at 2 a nd 6 cm. In a new detailed study, [Herbig & Grif fin (2006]l 
(see also Costero et al.[[2 006) discovered that 9^ Ori E is a double-lined spectroscopic 
binary consisting of two essentially identical mid-G-type components orbiting with a 
period of 9.89 days. This is a revision of the lower quality G-I-B5-B8 spectral type 
tabulated for component E in the Parenago] ([T954a ), Hillenbrand ( 1997), and Luhman 
et al. (20001 catalogs but originating from observations by |Herbig| ( il950| . Herbig & 
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Griffin also summarize clues that suggest significant optical variability of this system 
over the past 200 years. 



Ori F; Parenago 1892 Component F is a bright point source 4.5" SE (PA=120°) of 
9^ Ori C, appearing as a blue star in Figure 26 (see also 23 1; it was assigned a spectral 
type B8 by Herbig ( |1950 l. There are, however, no subsequent published observations 
to confirm the spectral type of this star. Existing high angular resolution studies imply 
this star is solitary down to a separation of 60 AU (Pe tr et al.|1998[ Simon et al.|1999[ l, 
which could perhaps explain why it is X-ray quiet ( |Stelzer et al.|2005' l. 



Regarding Components G, H During nineteenth century observations of the Trapez- 
ium, two additional sources within the Trapezium received designations as components 
of 6^ Ori . Components G and H are now known to be proplyds (See O'Dell et al.. Part 
II) corresponding to sources 2 and 3 of |Laques & Vidal, ( ,1979) . 



5. Variable Stars 



Photometric variability is a traditional technique for identifying young members of a 
known star forming region. Other methods such as kinematics (proper motions and ra- 
dial velocities), presence of lithium in the photosphere, and stellar/circumstellar activity 
exhibited as X-ray emission, UV and optical emission lines, and/or infrared excess have 
also been used. This section concentrates on variability as a young star selection tech- 
nique applied to Orion, in order to both provide historical perspective and connect to 
modem questions. 



Table 6. Modern variability surveys of the Orion Nebula. In addition to publica- 
tion details, the Table lists details of the observations such as whether photometry 
was published for the sources, the number of periodic stars recorded, and the ob- 
servatory used. Column entries with the value n.p. indicate that quantity was not 
published. 



Year 


Author 


Region 


Filter 


Phot? 




Np 


P(days) 


Obs. 


1990 


Walker 


ONC 


V 


Y 


5 


4 


0.4-3 


Lick 


1991 


Mandel 


Trap. 


Ic 


N 


150 


7 


6-14 


WO 


1992 


Attridge 


R < 0.25° 


Ic 


N 


525 


35 


2-17 


WO 


1995 


Eaton 


Trap. 


Ic 


N 


126 


11 


2-35 


WO 


1996 


Choi 


R < 0.25° 


Ic 


N 


525 


50 


2-20 


WO 


1999 
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Figure 27. Infrared light curves of variables in the Orion Nebula from Carpenter 
"etaL] ( |200T) . Left: star 2MASS J05342437-0452524 with a period of 8 days; Right: 
the eclipsing B star 9^ Ori Bi ^ (BM Ori); see also Sect.|4.4 
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Figure 28. Infrared Ught curves of irregular, long term variables in the Orion 
Nebula from ,Caipenter et al.| ( j2001^ . Left: Parenago 2171 (AO Ori); Right: Pare- 
nago 1617 (YY Ori). 



5.1. Optical & Infrared Variable Stars 

Historic Studies The Orion Nebula Cluster is a rich collection of variable stars. Ac- 



cording to Herbig ( 1982 1 the first variable star identified within the Orion Nebula was 



found by W.C. Bond in or around 1848, and is now known as AF Ori (Herbig 1982). 
Herbig notes in 1982 that in addition to the known eclipsing binary members of the 
Trapezium itself (9^ Ori A and 9^ Ori B) there were 17 named variables within the 
inner Orion Nebula Cluster (Trapezium region) and several hundred within the larger 
nebula. This area is known today to contain several thousand young stars, the vast ma- 
jority of which are known variables based on modern optical CCD or infrared array 
monitoring surveys with the cadence and sensitivity to detect them. 

The first extensive catalog of variable stars in (and near) the Nebula was derived 
(pretty much solely) by H. S. Leavitt using a large number (~ 20) of plates in the Har- 



vard collection. The catalog was published subsequently in Pickering & Leavitt ( 1904); 



it included approximately 70 variables and 30 suspected variables with minimum and 
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maximum brightness quoted. Additional notab le varia ble identification work was pe r- 
formed by |Brun| ( |l935| ), |Rosino| ( |l946[ |l956l ) and by jParenago & Kukarkin' ('1947'^ 
These and other stars were followed up spectroscopically, by e.g. Greenstein & Struve] 
( [T9431 ), [Herbigl ( fT950l ), and [Harol ( fT953] ). The analogy to "the objects of Joy" (1945), 
variable emission-line stars associated with bright or dark nebulosities and now known 
as T Tauri stars, was debated in these early works on the Orion variables, and eventually 
they were accepted as such, especially following photoelectric photometry studies such 
as that of |Walker| ( [19691 ). 

With improved technology, the variability could not only be identified and crudely 
characterized, but monitored in detail. [Walker ( 1990| ) undertook the first comprehensive 
study of rotational flux modulation in the Orion population, following up his earlier 
work on rotational broadening of spectral lines. He sought differences between the 
"weak" and "classical" emission stars as well as to understand the periods and nature 
of the surface disturbances causing the rotational modulation. He reported results on 5 
stars. 



Modem Monitoring Surveys With the advent of CCD detectors variability could be 
monitored and quantified for large numbers of members rather than just individual stars. 
In particular, there was interest in searching for periodic variables. Soon following 
the [Walker ( [1990[ ) paper was a series of papers by Herbst and collaborators (Mandel] 



& Herbst" 



99 1[ [Attridge & Herbst|[T992[ [Eaton et al.|[T995t [Choi & Herbst||1996' 



Table [6[ documents and compares the results from these modem survey beginning with 
Walker (1990). These works built up a database of periods in the Orion Nebula Cluster 
and presented the angular velocity distribution for the young stellar population. Of 
great interest was an apparent bimodal period distribution, which was interpreted in the 
context of "disk locking" whereby the slow rotators (with periods around 8 days) were 
assumed to be kept rotating slowly due to interaction between the stellar magnetosphere 
and the Keplerian accretion disk, while the more rapid rotators (with periods around 3 
days) were assumed free of such disk interactions. The gap or valley in the period 
distribution was interpreted as rapid evolution between the disk-regulated to disk-free 
scenarios. 

Additional wo rk on optical ([Stassun et all 1 999} [Herbst et al.[[20 00; R ebull|2001 



Herbst et al.||200T[ 12002; Stassun et al. 2006b, 20 07bUMariUi et al.||2007t jlrwin et al 



2007 1 and infrared (Carpen ter et al.,2001) variability broadened the discussion consid- 
erably. Included were not only the periodic stars with either cool spots (due to pho- 
tospheric inhomogeneities) or hot spots (due to accretion columns), but also eclipsing 
systems, and irregular variables dominated by accretion or extinction effects. |Herbst[ 



et al.[ ([2002[) claim that essentially every star (of 1500 monitored) is optically variable 

k-to-peak brightness variables at I- 
found the same mean peak-to-peak 



at the > 1 % level with half of the stars having peak-to-peak brightness variables at I- 
band larger than 0.2 mag. [Carpenter et al.[([2001 



amplitudes at near-infrared (JHK) wavelengths. 

The mass dependence and the activity/disk dependence of the periodicity was dis- 



cussed explicitly in several of the above works as well as in Stassun et al. (2001 2004a 



[2006b ) and RebuU et al. (2006 1. As a result of the young age and large amount of avail 



able rotation data, the ONC has become the de facto cluster setting the initial conditions 



Another significant study performed by P. P. Parenago was the development of the General Catalog of 
Suspected Variables with B. V. Kukarkin 
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for models of stellar angular momentum evolution. A compendium of the literature for 
the Orion Nebula and flanking field periodic stars is provided by 'RebuU et al. ( |2006[ ). 
Their results on variation in near-IR color with periodicity and long term irregularity 
provide clues to the origin of the variability, whether it be from time variable hot spots, 
from changes in line of sight extinction or due to changes in the geometry and rate 
of mass accretion. Near-IR light curves that illustrate different kinds of longer term 



variability are given in Figures 27 and 28 




2'10 4-10 6-10 B-IO 1-10 



12 
ID 



U 2 
D 



g A: 0.1 



1.5 



2.0 



Figure 29. Top: The Ughtcurve of the COUP 697 source, for which a modulation 
analysis yields an X-ray period similar to 10.2 days, as seen for the visible light 
starspots. Bottom: Here is the same lightcurve, af ter the flares are remove d, shown 
folded with the 10.2 day period. Reproduced from Flaccomio et al. (2005 1. 



5.2. X-ray Flares and Rotational Modulation of Variable Stars in Orion 

Concerning the geometry of the magnetic structures producing X-ray emission of T Tauri 
stars, solar-type coronal loops associated with multi-polar fields rooted in the stellar sur- 
face are probably the dominant source of the observed X-ray emission. The analyses 
of simultaneous optical and X-ray data for about 800 COUP stars show that optical and 
X-ray variability are very rarely time correlated ( Stassun et al.||2006b l, but the strong 



correlation between optical variability and X-ray luminosity is present (Stassun et al. 



2007b | ). This fits into the picture in which sights of optical variability may represent 
footprints of X-ray emitting coronal magnetic structures of complex magnetic topolo- 
gies. Studying 233 Orion-COUP PMS stars with known rotational periods, Flaccomio 



et al. ( |2005| ) detected X-ray rotational modulation in 16 stars, indicating that the stel 
lar surface has similar inhomogeneities in their photospheres and their coronae; it also 
suggests that the coronae in these cases are relatively compact R^,). An example 
of X-ray rotational modulation with the optical period is shown in Figure|29]for COUP 
source #697. In 7 other cases, Flaccomio et al.| (p005 ) find X-ray periods equal to half 
of the optical periods, suggesting two bright hemispheres in the X-ray corona. But 
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for a number of COUP stars with very powerful and hot (peak temperatures > 100 
MK) flares, Favata et al.'(2005) derive the length of the magnetic structures to be much 
greater than the stellar radius. These structures are probably too large to be stable, par- 
ticularly to centrifugal forces as the star rotates, in a solar-type geometry. They most 
likely extend from the star to the inner edge of a protoplanetary disk. 
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